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ABSTRACT 


This  report  is  concerned  with  the  theoretical  wave  resis¬ 
tance  of  an  air-cushion  vehicle  (ACV)  traveling  over  water  of 
uniform  finite  or  infinite  depth,  in  steady  or  unsteady  motion. 
Referring  first  to  steady  motion,  it  is  shown  that  the  un¬ 
realistic  oscillations  in  the  wave  resistance  curve  at  low 
Froude  numbers  found  by  previous  workers  can  be  eliminated 
by  using  a  smoothed  out  pressure  distribution  rather  than  one 
with  sharp  edges  studied  exclusively  in  the  past.  The  main 
result  of  unsteady  motion  calculations  is  that  the  peak  wave 
resistance  in  shallow  water,  even  in  moderately  accelerated 
motion,  is  appreciably  less  than  the  corresponding  steady- 
state  value.  In  fact,  cases  have  been  found  where  an  ACV 
starting  from  rest  under  the  action  of  a  constant  thrust  would 
seem  to  be  unable  to  cross  the  critical  depth  Froude  number 
on  the  basis  of  quasi-steady  estimates  of  wave  resistance, 
while  the  more  elaborate  unsteady  calculations  show  that  it 
has  sufficient  power  to  reach  its  final  supercritical  cruis¬ 
ing  speed.  An  interesting  feature  of  unsteady  motion  is  that 
besides  wave  resistance  there  is  another  mechanism  transfer¬ 
ring  energy  to  the  free  surface  which  is  here  called  the  dy¬ 
namic  sustention  power.  Contrary  to  intuition,  the  wave  re¬ 
sistance  in  unsteady  motion  over  finite  depth  sometimes  be¬ 
comes  negative  at  supercritical  Froude  numbers  before  finally 
approaching  zero  at  infinite  speed. 
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NOMENCLATURE 


half  length  of  craft 
1/2F 7 

frontal  area  of  craft 

half  width  of  craft 

velocity  of  ACV 

drag  coefficient 

water  depth 

propeller  diameter 

aerodynamic  drag 

momentum  drag 

Froude  number 

depth  Froude  number 

acceleration  due  to  gravity 

unit  vectors  in  x,  y,  z  direction 

propeller  advance  ratio 

wave  number  =  v//w2  +  u2 

fundamental  wave  number  =  g/c2 

propeller  torque  coefficient 

propeller  thrust  coefficient 

length  of  pressure  distribution 

Laplace  transform  operator 

mass  of  ACV 

mass  flow  into  air  cushion 


n 


=  unit  normal  to  free  surface,  directed  inwardly 
N  =  propeller  speed  of  revolutions 

p  -  cushion  pressure 

po  =  reference  cushion  pressure 

?  =  power  absorbed  by  propeller 

Pg,  P  ,  =  Kochin  functions  defined  by  Eq.  (3.14).  They  are 
Qg,  Qq  different  from  the  Kochin  functions  defined  by 

Wehausen  and  Laitone  (1960,  p.  558)  or  by  Eggers, 
Sharma  and  Ward  (1967)  ,  being  greater  by  the 
factor ,  cp  *cos  6  . 

P2D'^2D  =  Two"Dimensional  Kochin  functions,  cf.  Eq.  (3.28) 

Q  =  propeller  torque 

R  =  wave  resistance 

s  =  distance  traveled  by  ACV 

S  =  area  of  pressure  distribution 

t  =  time 

T  =  thrust  of  propulsor 

w,  u  =  induced  longitudinal  and  transv  se  wave  numbers 

W  =  weight  of  ACV  =  mg 

• 

WT  =  total  rate  of  work  done  on  water,  cf.  Eq.  (3.9) 

Wg  =  dynamic  sustention  power,  cf.  Eq.  (3.11) 

9 

WR  =  power  to  overcome  wave  resistance  =  cR 

x,  y,  z  -  coordinate  system  fixed  to  ACV,  cf.  Fig.  1 

* 

x  -  longitudinal  coordinate  in  fixed  reference  frame 
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e 
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=  longitudinal  cushion  fall-off  parameter 

-  transverse  cushion  fail-off  parameter 
=  /gk  •  tanh  (kd) 

=  water  shallowness  parameter 
=  free  surface  elevation 

-  propeller  efficiency 

=  wave  direction  with  respect  to  the  x  axis 
=  water  density 
=  air  density 
=  dummy  time  variable 

=  disturbance  velocity  potential  in  moving  frame 
such  that  v  =  V 4» 


=  n'th  term  of  asymptotic  expansion  for 
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SUPERSCRIPTS 
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=  stationary  frame  variable 
'  =  dummy  variable 
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"  at 

SUBSCRIPT 
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SUMMARY 


This  report  is  concerned  with  the  theoretical  wave  re¬ 
sistance  of  an  air-cushion  vehicle  (ACV)  traveling  over 
water  of  uniform  finite  or  infinite  depth,  in  steady  or  un¬ 
steady  rectilinear  motion. 

It  is  conventional  in  such  an  analysis  to  model  the 
ACV  with  a  given  pr<_  >*e  distribution  applied  to  the  free 
surface  of  an  inviscid  inconmressible  fluid  and  to  use  lin¬ 
earized  boundary  conditions  on  the  free  surface. 

The  results  obtained  by  this  approach  in  the  past, 
while  in  good  agreement  with  measurements  at  high  Froude 
numbers,  have  raised  two  questions  of  practical  signifi¬ 
cance.  First,  at  low  Froude  numbers  the  theory  predicts 
an  infinite  number  of  unrealistic  humps  and  hollows  in  the 
wave  resistance  curve.  Second,  when  the  depth  of  water  is 
small  compared  to  the  length  of  the  ACV,  the  steady-state 
peak  wave  resistance  at  the  critical  depth  Froude  number 
becomes  relatively  high  compared  to  the  wave  resistance  at 
the  cruising  speed,  which  is  typically  supercritical. 

It  is  shown  that  the  unrealistic  oscillations  at  low 
Froude  numbers  can  be  essentially  eliminated  by  using  a 
smoothed  out  pressure  distribution  in  contrast  to  the  sharp 
edged  distribution  used  by  previous  workers.  Moreover, 
for  a  rectangular  distribution,  this  effect  is  mainly  pro- 
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duced  by  the  smoothing  at  the  forward  and  after  edges. 

In  resolving  the  second  question,  the  primary  consider¬ 
ation  was  that  such  peak  resistance  can  only  represent  a 
transient  phase  in  the  practical  operation  of  any  ACV. 
This  naturally  suggests  tackling  the  unsteady  motion  prob¬ 
lem,  also  in  the  hope  that  in  accelerated  motion  the  ex¬ 
tremely  long  shallow  water  waves  may  never  have  enough  time 
to  build  up  to  their  peak  values. 

The  unsteady  theory  indeed  shows  that  at  reasonable 
accelerations,  an  ACV  can  pass  the  critical  depth  Froude 
number  without  encountering  unreasonably  high  wave  resis¬ 
tance  . 

It  is  seen  that  in  unsteady  motion,  besides  wave  re¬ 
sistance,  there  is  another  mechanism  transferring  energy  to 
the  free  surface.  This  is  the  dynamic  sustention  power 
and  represents  the  work  done  by  the  pressure  against  the 
relative  vertical  motion  of  the  free  surface  in  order  to 
maintain  the  altitude  of  the  ACV.  This  is  quite  indepen¬ 
dent  of  the  static  lift  power  required  to  support  the  air- 
cushion  vehicle  -  even  at  zero  speed. 

Results  of  several  sample  calculations  are  presented, 
including  many  for  a  two-dimensional  pressure  band  which 
relatively  easy  to  compute,  and  exhibits  the  phenomena  of 
interest  in  a  very  accentuated  manner.  Contrary  to  intu¬ 
ition,  in  two-dimensional  unsteady  motion,  the  wave  resis¬ 
tance,  and  even  the  total  rate  of  work  done  on  the  free 
surface  (including  sustention  power)  become  negative  at 
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some  supercritical  speed,  before  finally  approaching  zero 
at  infinite  speed.  However,  there  should  be  no  fundamental 
objection  to  this  phenomenon,  as  the  ACV  is  merely  recover¬ 
ing  some  of  the  energy  previously  expended  on  the  free  sur¬ 
face. 

It  is  found  that  the  Froude  number  at  which  the  maximum 
negative  wave  resistance  occurs  can  be  predicted  by  an  ap¬ 
plication  of  the  simpler  shallow  water  theory  -  in  which  the 
phenomenon  is  further  accentuated.  No  region  of  negative 
wave  resistance  was  encountered  for  a  three-dimensional 
pressure  distribution. 

Finally,  several  cases  of  the  inverse  problem  have  been 
calculated,  which  is  aimed  at  determining  the  velocity  pat¬ 
tern  for  an  ACV  starting  from  rest  under  the  action  of  a 
propulsor  of  given  thrust-speed  characteristics.  This  is 
treated  in  two  different  ways:  calculating  the  wave  resis¬ 
tance  in  a  truly  unsteady  manner,  and  on  the  simplified 
quasi-steady  basis.  All  other  components  of  drag  are  as¬ 
sumed  to  be  strictly  quasi-steady.  The  results  show  that 
the  shape  of  the  propeller  thrust  and  torque  coefficient 
curves  has  only  a  minor  effect  on  the  velocity  pattern.  On 
the  other  hand,  the  effect  of  overloading  the  ACV  is  found 
to  have  crucial  effects  on  its  ability  to  surpass  the  criti¬ 
cal  depth  hump. 

In  this  respect,  the  simpler  quasi-steady  calculations 
lead  to  unnecessarily  pessimistic  estimates  of  the  velocity 
pattern.  Under  certain  circumstances  in  relatively  shallow 
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water,  the  quasi-steady  analysis  would  suggest  that  the  ACV 
could  not  overcome  the  critical  hump  with  the  available  pow¬ 
er,  while  the  more  elaborate  unsteady  calculations  show  that 
it  has  indeed  adequate  power  to  reach  its  final  cruising 
speed. 
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1  -  INTRODUCTION 


1.1  -  PREVIOUS  WORK:- 

The  hydrodynamic  aspects  of  an  air-cushion  vehicle  (ACV) 
can  be  studied  by  assuming  its  action  to  be  equivalent  to  that 
of  a  pressure  distribution  acting  on  the  free  surface  of  the 
water.  This  idealization  neglects  any  physical  contact  of 
the  lowar  edge  of  the  craft  with  the  water.  It  also  assumes 
that  the  flow  of  air  escaping  under  the  periphery  is  inviscid, 
and  therefore  produces  no  spray. 

Havelock,  in  some  of  his  early  papers  (1909  ,  1914  $.  1926) 
was  the  first  to  treat  the  theoretical  problem  of  the  wave 
resistance  of  a  pressure  distribution.  His  interest  in  pres¬ 
sure  disturbances  lay  in  a  desire  to  represent  the  motion  of 
a  ship.  As  a  results,  most  of  the  distributions  that  he  chose 
to  analyse  were  very  smooth  and  were  not  typical  of  the  ACV. 

However,  later  on,  Havelock  (1932)  derived  the  general 
expression  for  a  pressure  distribution  traveling  at  a  constant 
speed.  In  this  paper,  he  also  found  the  relationship  between 
the  pressure  acting  on  the  free  surface,  and  the  equivalent 
source  distribution. 

Lunde  (1951a)  extended  the  theoretical  treatment  to  cover 
the  case  of  an  arbitrary  distribution  moving  over  finite  depth. 
Other  workers  have  obtained  numerical  results  for  pressure 
distributions  which  are  directly  applicable  to  the  ACV.  These 
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include  Newman  and  Poole  (1962)  who  considered  the  case  of 
motion  in  a  restricted  waterway  such  as  a  canal.  They  cal¬ 
culated  the  two  cases  of  a  constant  pressure  acting  over  a 
rectangular  area,  and  over  an  elliptical  area.  The  most 
striking  feature  of  their  results  is  the  very  strong  inter¬ 
action  between  the  bow  and  stern  portions  of  the  distribut¬ 
ion.  Particularly  for  the  rectangular  distribution  (where 
the  interaction  would  be  greater) ,  there  are  displayed  a 
series  of  humps  and  hollows  in  the  resistance  curve.  A 
hump  occurs  when  the  bow  and  stern  wave  systems  are  in  phase 
and  combine  to  give  a  trailing  wave  of  a  maximum  height.  A 
hollow  occurs  when  the  two  wave  systems  are  out  of  phase  by 
half  a  wavelength  giving  a  combined  amplitude  of  a  minimum 
height. 

The  interference  effects  are  found  to  be  stronger  for 
large  beam  to  length  ratios,  as  would  be  expected  from  the 
above  argument,  since  the  wave  motion  becomes  more  nearly 
two-dimensional  for  a  wide  craft. 

Barratt  (1965)  also  computed  the  wave  resistance  of  rec¬ 
tangular  and  elliptical  pressure  distributions,  but  for  the 
case  of  unrestricted  water.  His  results  are  to  some  extent 
similar  to  those  for  the  canal.  In  deep  water,  the  main  (or 
"last")  hump  occurs  at  a  Froude  number  given  by  F  =  1//F. 

In  water  of  finite  depth  this  hump  is  shifted  to  a  lower 
Froude  number,  and  for  sufficiently  shallow  water  occurs  at  a 
depth  Froude  number,  F^,  equal  to  unity  (i.e.,  at  the  critical 
speed).  One  difference  between  these  two  sets  of  results  is 
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pointed  out  by  Hewman  and  Poole.  For  a  canal  of  finite 
width,  the  theory  predicts  a  discontinuity  in  the  wave  re¬ 
sistance  at  the  critical  depth  Froude  number.  The  resistance 
is  higher  just  below  the  critical  speed  than  just  above  it. 
However,  for  an  infinitely  wide  canal,  there  is  no  discon¬ 
tinuity,  but  there  is  a  sudden  change  in  slope  at  the  critical 
speed . 

Havelock  (1922)  also  presented  some  results  for  a  very 
smooth  pressure  distribution  over  water  of  finite  depth. 

These,  too,  clearly  show  the  shift  of  the  main  hump  and  the 
increase  in  its  magr itude  in  shallower  water.  Havelock's 
curves  display  only  the  main  hump.  The  secondary  and  other 
humps  do  not  occur  because  of  his  choice  of  pressure  distri¬ 
bution  . 

Recently  a  number  of  experimental  programs  have  been 
carried  out  in  order  to  check  the  above-mentioned  theoretical 
results.  Chief  workers  in  this  field  are  Everest  (1966a, 

1966b  &  1967)  and  Hogben  (1966a).  The  main  question  pointed 
out  in  these  papers  is  the  resolution  of  the  total  drag  on  the 
ACV  into  its  components.  These  components  are  often  consider¬ 
ed  to  be:  Wave  Resistance,  Aerodynamic  Drag,  Momentum  Drag 
dnd  Water  Contact  Drag. 

The  aerc dynamic  drag  (or  profile  drag)  is  assumed  to  be 
that  resistance  acting  on  the  model  if  it  were  tested  in  a 
wind  tunnel  with  the  engines  not  running. 

The  momentum  drag  is  that  due  to  the  change  in  direction 
of  the  air  supplying  the  cushion  as  it  enters  the  fan  intakes 
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In  fact,  we  should  consider  two  components  here:  Inlet  Mo¬ 
mentum  Drag,  and  Outlet  Momentum  Drag.  The  outlet  moment¬ 
um  drag  is  associated  wi th  the  changes  in  direction  and  ve¬ 
locity  of  the  air  as  it  escapes  from  the  cushion,  and  can  be 
either  positive  or  negative,  depending  on  the  trim  of  the 
craft. 

The  water  contact  drag  is  due  to  any  touching  of  the 
lower  edge  of  the  craft,  or  of  the  skirts,  with  the  water. 

Due  to  the  extremely  non-linear  nature  of  this  effect,  the 
drag  (or  possibly  thrust)  due  to  spray  from  the  cushion  hit¬ 
ting  the  craft  is  usually  included  with  it. 

While  the  first  thr^e  resistance  components  defined 
above  may  be  studied  from  a  theoretical  approach,  the  water 
contact  drag  only  lends  itself  to  an  experimental  study.  To 
this  end,  Everest  (1966a)  estimated  the  water-wetting  resis¬ 
tance  by  eliminating  it  -  using  a  thin  polythene  sheet  floa¬ 
ting  on  the  water  surface.  This  technique,  however,  intro¬ 
duces  the  question  about  the  tensile  forces  in  the  sheet. 

The  resistance  breakdown  is  further  discussed  by  Hogben  (1966a) 
where  he  provides  a  careful  definition  of  each  component. 

These  experimental  results  were  generally  obtained  by  a 
dynamometer  measuring  the  total  drag.  Then  the  aerodynamic, 
momentum  and  water  contact  drag  components  were  estimated 
and  subtracted  in  order  to  make  a  comparison  with  the  theore¬ 
tical  wave  resistance.  The  agreement  appears  to  be  quite 
good  in  regard  to  the  range  in  values  of  the  wave  resistance. 
However,  the  large  scatter  in  the  data  makes  it  difficult  to 
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draw  precise  conclusions.  The  authors  suggest  that,  in 
addition  to  the  main  hump,  they  can  detect  a  secondary  and 
possibly  a  third  one,  and  that  these  are  out  of  phase  with 
the  theoretical  humps  by  no  more  than  0.05  on  the  Froude  num¬ 
ber  scale. 

An  explanation  given  for  the  non-appearance  of  more  humps 
is  based  on  the  fact  that  the  lower  speed  humps  predicted  by 
the  linearized  theory  correspond  to  a  wave  pattern  whose  max¬ 
imum  slope  is  too  large  from  physical  considerations.  Hogben 
gave  a  two-dimensional  argument  (1965)  showing  that  the  max¬ 
imum  ratio  of  wave  height  to  length  is  about  1/7.  This  would 
preclude  the  development  of  any  humps  above  the  third  or 
fourth  (depending  on  the  cushion  pressure) . 

Further  experimental  work  by  Everest,  Willis  and  Hogben 
(1968  &  1969)  dealt  with  the  wave  resistance  of  an  ACV  at  an 
arbitrary  angle  of  yaw.  This  problem  was  also  studied  theore¬ 
tically  by  Mui chy  (1970).  In  these  experiments,  the  wave 
resistance  was  measured  directly  from  the  wave  pattern.  As 
a  result  there  is  less  scatter  in  the  data  since  the  rather 
doubtful  technique  of  estimating  the  wetting  drag  is  elimin¬ 
ated.  The  experimental  results  here  are  generally  low  com¬ 
pared  with  the  theory,  the  difference  being  usually  limited 
to  about  10%  but  is  occasionally  as  much  as  50%  at  certain 
Froude  numbers. 

The  outcome  of  these  investigations  is  that  the  main 
hump  drag  is  relatively  large  in  relation  to  the  installed 
propulsive  power  of  typical  ACVs.  In  addition,  the  prob- 
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lem  is  more  acute  ir.  shallow  water  where  the  thrust  margin  of 
some  craft  has  been  found  insufficient  to  surpass  the  hump. 

It  has  been  found  necessary  from  speed  and  economic  considera¬ 
tions  to  operate  in  the  cruising  condition  at  a  Froude  number 
of  at  least  1.3.  This  is  well  above  the  hump  speed  and  the 
wave  resistance  is  accordingly  smaller. 

However,  in  a  real  situation,  thr  craft  does  not  operate 
steadily  at  the  hump  speed.  In  fact,  the  procedure  is  to  ac¬ 
celerate  through  it  as  quickly  as  possible.  Under  a  non¬ 
steady  condition  it  appears  quite  reasonable  to  anticipate 
that  the  large  amplitude  wave  pattern  at  the  hump  speed  will 
not  have  time  to  establish  itself  -  thus  leading  to  a  less 
aggravated  problem  in  shallow  water.  Some  experiments  on  a 
rectangular  model  by  Everest  (1966b)  confirm  this.  Under 
certain  conditions  in  finite  depth  he  has  found  a  reduced 
resistance  peak. 

1.2  -  PRESENT  STUDY : - 

These  considerations  point  towards  a  theoretical  inves¬ 
tigation  of  the  wave  resistance  during  accelerated  motion. 
Already,  the  problem  for  a  ship  has  been  treated  by  Sretensky 
(1939),  I  unde  (1951b,  1953a  &  1953b)  and  Shebalov  (1966). 

These  workers  have  derived  the  linearized  result  for  the  re¬ 
sistance,  but  produced  no  computed  values.  Wehausen  (1964) 
computed  the  resistance  of  a  ship  model  with  a  constant  ac¬ 
celeration  from  rest  up  to  a  speed  which  was  then  held  fixed. 


7 


His  results,  however,  consist  of  asymptotic  expressions  for 
large  values  of  the  time,  and  there  are  no  data  for  the  re¬ 
sistance  during  the  acceleration  phase  of  the  motion.  Wehau- 
sen's  interest  stemmed  from  a  desire  to  know  the  required 
length  of  a  model  ship  test  before  the  steady  state  resistance 
is  achieved. 

With  regard  to  the  unsteady  motion  of  a  surface  pres¬ 
sure  distribution,  Havelock  (1916)  has  computed  the  resistance 
of  two  particular  distributions  whose  motion  is  suddenly  es¬ 
tablished  from  rest,  and  then  continued  at  a  constant  speed. 

On  the  other  hand,  Djachenko  (1966)  has  derived  an  expression 
for  the  resistance  of  an  arbitrary  pressure  distribution  for 
a  general  acceleration  pattern.  He  also  presented  some  re¬ 
sults  for  a  two-dimensional  distribution. 

The  transient  problem  is  of  interest  mainly  for  two  rea¬ 
sons.  First,  since  a  numerical  solution  to  this  problem  will 
make  it  possible  to  calculate  the  resistance  history  as  a  func¬ 
tion  of  the  acceleration  pattern,  one  could  expect  to  find  op¬ 
timum  acceleration  programs  that  would  reduce  the  peak  wave 
resistance,  or  the  peak  total  resistance. 

Second,  it  is  of  interest  to  examine  the  inverse  prob¬ 
lem:  the  resulting  motion  under  the  action  of  a  propulsive 

device  of  given  thrust-speed  characteristics.  This  is,  of 
course,  the  more  natural  problem  and  any  theoretical  results 
could  also  be  checked  by  experiment. 

It  is  believed  that  there  is  more  justification  for  ap¬ 
plying  the  linearized  theory  to  the  transient,  rather  than 
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the  steady,  problem  because  of  the  diminished  likelihood  of 
an  excessively  steep  wave  system  building  up. 


2  "  THE  POTENTIAL  FUNCTION 


2.1  -  PROBLEM  STATEMENT:- 

The  air-cushion  vehicle  will  be  represented  by  a  pres¬ 
sure  distribution  p(x,y)  acting  on  the  free  surface,  and 
traveling  with  the  speed  of  the  craft.  Two  right-handed 
coordinate  systems  will  be  used,  as  shown  in  Fig.  1.  The 
system  xyz  moves  with  the  craft,  z  being  vertically  up¬ 
wards  and  x  being  in  the  direction  of  the  rectilinear  mo- 

* 

tion.  The  second  axis  system  x  yz  is  fixed  in  space. 

The  relationship  between  the  coordinates  is  then  given  by 

* 

X  =  x  -  s  (t) 


* 


X 


t 

■ 

c(t)  di, 

0 


(2.1) 


where  c  and  s  are  the  velocity  of  the  model,  and  its  dis¬ 
tance  traveled,  respectively. 

The  velocity  potential  satisfies  the  Laplace  equation: 

*  2  * 

v  <t>  -  0 ,  (2.2) 

* 

where  <J>  is  the  velocity  potential  in  the  stationary  frame 

(such  that  the  velocity  is  the  positive  gradient  of  the  po- 

* 

tential) ,  and  V  also  implies  differentiation  in  this  frame. 
If  4)  is  the  perturbation  potential  in  the  moving  frame,  we 


9 


10 


may  say  that 

it  * 

$ (x,y ,z,t)  ~  0  (x  ? y / z / 1 ) i  (2.3) 

* 

if  x  and  x  correspond  to  the  same  point  in  space.  From 
Eq.  (2.1)  we  have 


3 


3x 


3_ 
3x  ' 


(2.4) 


so  that  Eq.  (2.2)  may  be  written  as 


V2$  =  0  . 


(2.5) 


The  kinematic  boundary  condition  on  the  free  surface  re¬ 
quires  that  a  particle  on  it  remains  there.  That  is. 


D 

Dt 


C (x,y,t) 


u 


o , 


(2.6) 


where  £  is  the  elevation  of  the  surface.  In  terms  of  the 
moving  coordinates  we  have 


D 

Dt 


"■k +  ,u 


-  * 
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where  u,  v  and  w  are  the  perturbation  velocities.  Com- 
bininc  this  with  the  exact  kinematic  condition,  Eq.  (2.6), 
and  substituting  fcr  u,  v  and  w: 

-  Ocx-*y?y]  z.c  *  °‘ 


At  this  stage  the  second  order  terms  may  be  dropped,  and  the 
remaining  terms  written  as  a  Taylor  expansion  about  the 
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point  z  =  0: 


5  + 


C.  +  c  c 
t  x 


0  . 


Finally,  the  linearized  kinematic  condition  on  the  free  sur¬ 
face  may  be  obtained  by  dropping  the  higher  order  quantities 
again: 


?t  +  c 


=  0. 


(2.7) 


The  dynamic  condition  on  the  surface  -  the  Bernoulli  equa¬ 
tion  -  in  terms  of  the  stationary  coordinates  is 


[♦*t 


+ 


iu*  2 

x* 


+  <J> 


+  4>  ,  > 


u 


* 

£  + 
p 


95  =  f,  (2.8)) 


where  p  is  the  water  density  and  g  is  the  acceleration  due 
to  gravity,  while  f  is  an  arbitrary  function  of  time. 

Again,  we  may  drop  the  squared  terms  as  these  are  of  higher 
order.  In  addition,  f  may  be  put  identically  to  zero  with¬ 
out  loss  of  generality.  Also,  from  Eqs.  (2.1)  and  (2.3)  we 


<P 


* 

t 


c  <p 


have 
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Thus  Eq.  (2.8)  reduces  to 

[*t  '  0  *x]z=c  = 

As  before  we  nay  expand  the  first  term  in  a  Taylor  series 
about  z  =■  0  ,  and  then  drop  the  higher  order  terms  to  ob¬ 
tain  the  linearized  dynamic  condition: 

[>  •  C  ♦xJz-O  +  p  +  «  “  0  •  <2-9) 

The  combined  free  surface  condition  is  obtained  from 
Eqs.  {2.7)  and  {2.9)  after  eliminating  c  : 

[*tt  +  °'*xx  -  2c  +xt  -  ®  »x  +  9  *z]z=0  "  c  Px/o  •  (2-10> 

The  last  boundary  condition  to  be  satisfied  is  that 
there  is  no  flow  thri  ugh  the  water  bed: 

[*z]z— d  =  0  <  f2-11’ 

where  d  is  the  depth  of  the  water. 
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2.?  -  SOLUTION  FOR  THE  POTENTIAL;- 

The  solution  of  this  set  of  equations  can  be  obtained  by 
an  application  of  the  double  Fourier  transform  pair: 


oo  co 


f(w,u)  = 


dx 


dy  f(x,y)  exp{-i(wx  +  uy)  ) 


_co  —  oo 


(2.12) 


CD  OO 


and  f  { x ,  y )  =  ~ 


dw 


du  f(w,u)  exp  (i(wx  +  uy)  )  , 


•00  — OO 


and  the  Laplace  transform  pair: 


f  (q)  =  f  (t)  exp(-qt)  dt 
0 

6+i°° 

• 

and  f  (t)  =  f  (q)  exp(qt)  aq,  (2.13) 

5“i“ 

6  being  a  positive  constant. 

Using  the  rules  for  transforming  a  derivative,  the  La¬ 
place  equation  (2.5),  under  the  transformation  (2.12),  be¬ 


comes  : 
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(iw)  2 <f>  +  (iu)  2<J>  +  <t>zz  =  0  , 


where 


<j>  =  4>(w,u;z,t) 


Hence 


*zz  -  k‘*  ~  0  ' 


where 


k2  =  w2  +  u2  . 


(2.14) 


A  solution  for  <J>  is 


<J>  =  A(w,u;t)  cosh(kz  +  B(w,u;t))  ,  (2.15) 


where  A  and  B  are,  at  the  moment,  arbitrary.  B  may  be 
found  by  transforming  the  bed  condition,  Eq.  (2.11),  giving 


[♦  J«-d  ■ 0  - 


and  substituting  Eq.  (2.15).  Thus 


A*k* sinh  (kz  +  B) 


]z=-d  ‘  0  ' 


giving 


B  =  kd 


(j>  =  A  cosh  (k  (z  +  d) ) 


(2.16) 
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Applying  now  the  Fourier  tvansfoon  to  the  combined  free 
surface  condition,  Eq.  (2.10),  and  using  $  given  by  Eq. 
(2.16),  we  obtain 


|^tt  •  cosh  (k  (z  +  d)  )  +  c2  (iw)  zA*cosh  (k  (z  +  d)  ) 


-  2c* iw* A^1 cosh (k (z  +  d))  -  c iw*A’cosh (k (z  +  d) ) 
+  gk*sinh(k(z  +  d))l  z_q  =  ciwp/p  . 


Or,  more  simply, 

Afct  -  2icw*Afc  +  A ( y 2  -  c2w2  -  icw)  =  icw*  sech  vkd)  •  ,  (2.17) 

in  which  y2  -  gk  tanh(kd)  .  (2.18) 

A  substitution  that  makes  the  coefficients  in  Eq .  (2.17) 
constant  has  been  found  by  Lunde  (1951b) .  It  is 

*v  «•# 

A  (w,u;  t)  =  y(w  exp(iws(t))  .  (2.19) 

After  putting  this  in  Eq.  (2.17),  and  some  simplification, 
one  obtains 

y  +  y2y  =  icw*  sech  (kd) -2.»exp  (-iw- s  (t)  )  . 
tt  P 


16 


We  now  take  the  Laplace  transform,  Eq.  (2.13),  of  the 
previous  equation  to  get 

(qz  +  y2)  X  =  iw  *3ech  (kd)  •£' z,  (c  *exp(-iw*s  (t) )  )  , 

where  L  is  the  Laplace  transform  operator.  Thus 

X  =  i  •^■•sech  (kd)  •  L(sin  yt )  * L  (c  *exp  {-iw  •  s  (t)  ))  . 

Using  the  convolution  theorem  to  perform  the  inverse 
Laplace  transform  of  this,  we  obtain 


X  = 


i  .iiR.  sech  (kd) 
PY 


c ( i ) • sin (y  (t 


t) ) *exp  (-iw*s (i) )  di  . 


Eqs.  (2.19)  and  (2.16)  are  now  used,  together  with  the 
above  equation,  to  yield 


<(>  = 


wp  cosh (k (z  +  d) ) 
1* py*  cosh (kd) 


c (x ) • sin (y  (t  -  t ) ) 


•exp(iw(s(t)  -  s  (  t  ) )  '  di  .  (2.20) 

The  double  Fourier  transform,  Eq.  (2.12),  can  be  used 
to  express  p  in  Eq.  (2.20)  in  terms  of  p,  and  then  the 
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inverse  transform  is  taken  to  give  the  disturbance  potential, 
<1?  ,  in  terms  of  the  moving  coordinate  system,  xyz  : 


Here  k2  =  w2  +  u2 

as  before. 

and 

P' 

=  p1  (x' ,y 1 )  ,  defined 

over  the  area  S'  , 

while  x' 

and 

y' 

are  dummy  coordinates 

in  the  moving  reference  frame. 


3  -  THE  WAVE  RESISTANCE  AND  SUSTENTION  POWER 


3,1  -  DEFINITIONS 

In  this  section  expressions  for  the  wave  resistance  will 
be  derived.  Also  the  computer  studies  showed  that  under  cer¬ 
tain  circumstances  this  resistance  could  become  negative. 

Hence  it  is  of  interest  to  examine  the  total  rate  of  energy 
input  to  the  water,  as  well. 

Using  the  stationary  frame  of  reference,  the  wave  resis¬ 
tance  may  be  defined  as  the  horizontal  component  of  the  cush¬ 
ion  pressure  force  acting  on  the  free  surface.  Thus 


R 


*  *  *  * 
p  (x  ,y,t)  C  x*  dx  dy  , 

* 


★ 

where  the  superscript  indicates  that  the  variable  is 

given  in  terms  of  the  fixed  frame.  This  formula  may  be  ex¬ 
pressed  in  the  moving  frame  as 
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Now  the  power  input  to  the  water  (in  the  fixed  frame) 
is 


W„ 


*  *  *  *  * 
p  (x  ,y,t)  n  *v  ds  , 

.  i 
* 

S 


(3.2) 


where  n  is  the  unit  normal  vector  at  the  free  surface  and 
directed  inwardly,  and  v  is  the  velocity  of  the  particles 
of  water  on  the  surface.  By  definition 


5*  =  U*x*  i  +  C*y  j  "  k)/A  +  U*x*)  2  +  U*y)  2  (3.3) 

and  v  ~4’x*i+(f,yj  +  (},zk,  (3.4) 

In  the  stationary  reference  frame,  the  kinematic  con¬ 
dition  on  the  free  surface  is 


D 

Dt 


(z 


C* (x* ,y , t) ) 


* 


0  , 


or 


[' 


-  G 


* 

4>  * 


~  C 


-  C 


tl  z=i 


*  =  0 


(3.5) 


In  addition,  we  have 
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*  /" 

.  J  _  _  /*! 


dS  =  dx  dy  /l  +  U  *)  ‘  +  U  2  . 

x  y 


(3.6) 


If  we  now  combine  Eqs.  (3.2)  to  (3.6),  we  obtain 


wT  =  - 


*  *  *  * 

P  (x  /Y/t)  C  t  dx  dy 


(3.7) 


A  transformation  to  the  moving  frame  is  now  in  order, 
and  using  the  relation 


* 

C 


t 


’x 


(3.8) 


the  final  expression  for  the  total  rate  of  work  done  on  the 
water  results: 


Note  that  the  expressions  for  the  resistance  and  total 
power  (Eqs.  (3.1)  and  (3.9))  are  exact  -  they  have  not  been 
linearized.  The  total  power  contains  the  power  required  to 
overcome  the  wave  resistance: 
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W  =  c  R  .  (3.10) 

K 

The  other  part  of  the  total  power  we  shall  call  the  sus¬ 
tention  power: 


W„  =  - 


p(x,y)  ;  dx  dv 


(3.11) 


Thus  the  sustention  power  represents  the  rate  at  which 
work  is  done  by  the  cushion  pressure  against  any  vertical 
motion  of  the  free  surface,  and  is  positive  for  an  average 
downward  movement.  This  power  term  should  not  be  confused 
with  the  power  required  to  provide  the  air  cushion.  The 
sustention  power  would  come  from  changes  in  gravitational 
potential  energy  of  the  ACV,  as  it  heaves  and  trims  during 
the  unsteady  motion.  During  steady  motion,  the  position  of 
the  free  surface  does  not  move  relative  to  the  craft.  In 
this  case  the  sustention  power  is  zero. 
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3.2  -  WAVE  RESISTANCE: - 

We  use  Eq.  (3.1)  defining  the  resistance,  and  the  dy¬ 
namic  cor.dition  on  the  free  surface,  Eq.  (2.9). 


R 


P 


z=0 


-  Px/e 


dx  dy  . 


After  integrating  with  respect  to  x  ,  the  last  term  in 
the  integrand  gives 


which  is  zero  for  any  typical  distribution.  We  may  now  sub¬ 
stitute  the  expression  for  4>  given  by  Eq.  (2.21),  and  af¬ 
ter  some  reduction,  obtain 


' 

■ 

-  c 

0  CO 

’ 

47T2pg 

p  dS 

P' 

dS’ 

C(T)  dT 

dw 

c 

;  £ 

•  1 

1 

•  •  4 

0  “CO  —oo 

w2  •  cosp gk*  tanh  (kd)  •  (t 

-  T  ) 

• 

•  exp^i  (w  (x 

:  -  x'  + 

s  (t) 

-  s  (x 

) )  +  u  (y  ■ 

-  y  * )  > 

du 


(3.12) 
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Since  we  only  require  the  real  part  of  this  equation,  we 
rewrite  the  exponent  factor  as 

factor  =  cosf"w(x  -  x’  +  s(t)  -  s(x))  +  u(y  -  y')l 


=  cos 

w  (x  - 

x '  ) 

+  u  (y  - 

■y,] 

•cos 

w(s  (t) 

-  s  {  t)  )  J 

-  sin 

w{x  - 

X') 

+  u  (y  ■ 

-  y') 

•sin 

w  (s  (t) 

-  s (t) ) J 

The  second  term  on  the  right  Land  side  is  odd  with  respect 
to  the  pair  x  -  x'  and  y  -  y'  ,  and  therefore  contributes 
nothing  to  the  integral.  Expanding  the  first  term,  now, 


factor 


- 


ros(w(x  -  x'))*cos(u{y  -  y')) 


-  sin(w(x  -  x'))*sin{u(y  -  y 1 )  )  J*  cos^w  (s  (t)  -  s(t))J 

This  time  the  second  term  is  odd  with  respect  to  w  and  u  , 
and  therefore  gives  no  contribution.  We  may  expand  the  fac¬ 
tor  a  third  time  to  yield 


factor  =  |^cos  <wx)  *cos  (wx  ' )  +  sin  (wx)  •  sin  (wx  1 )  > 


{continued  over) 
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'^cos  (uy)  *cos  (uy ' )  +  sin  (uy)  *sin  (uy  ’ )  Jx 


>  w 


>cos  w  (s  (t)  -  s  (  t)  ) 


]• 


If  we  place  this  factor  back  into  Eq.  (3.12),  the  de¬ 
sired  expression  for  the  resistance  is  obtained: 


IL  is  interesting  to  compare  Eq.  (3.13)  with  the  ex- 


pression  obtained  by  Lunde  (1951b)  for  a  ship.  An  addi¬ 
tional  term  consisting  of  a  fourfold  integral,  and  being  sim¬ 
ply  proportional  to  the  instantaneous  acceleration,  occurs 
in  his  formula.  It  represents  a  type  of  added  mass,  and  is 
zero  if  the  singularity  distribution  lies  on  the  free  sur¬ 
face.  This  is  the  case  of  a  pressure  distribution. 


3.3  -  TOTAL  POWER 

The  total  power  is  defined  by  Eq.  (3.9).  Using  the 
kinematic  condition,  Eq.  (2.7),  the  power  may  be  written  as 


wT  = 


|  ?  [*2]z-0  dx  dy  • 


The  equation  for  <f>  ,  Eq.  (2.21),  is  now  substituted  to 


give 


00  00 


wm  = 


-1 


T  4r2pg 


p  ds 


p'  as’ 


.  j 

S' 


c(x)  di 


dw 


du 


“CO  “  oo 


w* /gk*  tanh (kd) *  sin  •tanh(kd) • (t  -  t) J< 


-jterrr:*-?-  -  »'  ■  1  iiy i>»  t -^- --' '-  -  - 


2  6 

'I’i.i.  '.‘S|jri.‘!'.ii  .ion  may  now  be  simplified  in  a  manner  iden¬ 
tical  to  that  used  on  ttq .  (3.12)  for  the  resistance.  We 

then  obtatn  tor  the  total  rate  of  work  done  on  the  water  by 

tin-  pres. -nie: 


•  r 

1  :  ...  | 
"T  S-7u  :  ^ 

r 

dw  |  du  w  •  I'gk  •  tanh  (kd)  x 

1 

l  ! 

i>  0 

6 

-j  p  ■’  i  '  ■'  -t  o  •’  +  c 

i  < >  t :  t  * 

v]- 

“sin  v  g k • tank  (kd)  •  ( t 
L 

-  i)  -sin  w(sft)  -  s(t)  )  .  (3.15) 

3.4  -  1 '1  y i  i A VAT A  WAVE 

RESISTANCE  :  - 

The  steady  state  resistance  may  be  obtained  from 


F.q.  i3  .  1 J  )  ey  letting  the  velocity  of  the  model  be  constant 

for  a  '  os' ,  •  :  r' ex-amp  1  c ,  if  we  take  the  velocity  of 

tin  ■  < '  i:  i:  •  ‘  .  •  suddenly  os  tab]  Lshed,  and  then  fixed  at  a 

value  •  ,  !  :  >•  t  ir.K-  i  ntogra]  in  Kq .  (3.13)  may  be  carried 

O'.i  t  ana  !  t  >  • .  '  to  ;;  i  vo 


BEST  AVAILABLE  COPY 
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R  = 


2  7T  2  pg 


w2  dw 


du  p  2+p2+Q2+Q2  * 
e  *0  e  o 


sin(y  +  wc)t 
y  +  wc 


sin(y  -  wc)t 
y  -  wc 


(3.16) 


with  y2  =  gk«tanh(kd) 


=  g  /u  2  +  w  2  •  tanh  ( /u  ^  +  w2*d)  , 


as  before. 

We  want  to  know  the  value  of  R  for  very  large  time. 
As  t  -+  oo  ,  the  oscillations  (about  zevo)  due  to  the  sine 
terms  increase,  so  thct  there  is  no  contribution  to  the  in¬ 
tegral  except  when 


y  -  w  c  =  0  .  (3.17) 

Eq.  (3.17)  gives  the  relationship  between  the  trans¬ 
verse  and  longitudinal  wave  numbers  for  waves  which  travel 
at  the  speed  of  the  model.  We  may  now  quote  the  following 
result  from  Wehausen  (1964)  or  Wehausen  and  Lai tone  (I960, 
p.  477): 


sin  ( :<  -  x  Jt 


f  (x  o )  if  a  <  x  o  <  b 


It  is  required  here  that  f  is  smooth  in  the  neighborhood 
of  x  -  x o  ,  and  that  f(x)/x  is  absolutely  integrable  in 
the  range.  The  theorem  below  may  then  be  derived  from  this 
result : 


f(x)  5in  dx 

'  '  g  (x)  ax 


=  tt  l  f  (xi)/|g1  (Xi)  |  ,(3.18) 


where  are  the  zeros  of  g(x)  ,  assuming  that  a  <  x^  <  b 

and  g'  (x^)  f  0  . 

First  we  rewrite  Eq.  (3.16)  in  terms  of  polar  coor¬ 
dinates  : 


w  =  k  cos  0  , 


u  =  k  sin  0  , 


(3.19) 


^  (w,u)  _ 

TTkToT  "  k  ' 


and  then  apply  the  second  theorem  to  the  k  integral.  It 


is  seen  that  if 


g  (k ,  ■•■) 


w  c 
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or 


g(k,0)  =  /qk‘  tanh(kd)  -  ck*ccs  8  , 


then 


g (k, e) 


=  2°' 


•cos  0  1  -  kod*sec20-sech2 (kd) 

.. 


]• 


where  k  is  the  solution  of 


k  -  k0  sec2 0 • tanh (kd)  =  0  ,  (3.20) 


with 


k0  =  g/c2  . 


(3.21) 


If  we  now  make  the  necessary  substitutions  detailed 
above,  Eq.  (3.16)  reduces  to  the  following  single  integral 
for  the  resistance  of  an  \CV  during  steady  motion: 


Here  the  lower  limit  for  0  is  taken  as  0i  ,  the 
smallest  positive  value  of  0  which  can  satisfy  Eq.  (3.20) 
for  a  real  k  .  It  is  given  by 
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0  for  k  0d  ->  1 

=  .irccos  .  1c  oil  for  k0d  <  1 


(subcritical  speed) 

(3.23) 

(supercritical  speed) 


Eq.  (3.22)  is  the  same  as  that  given  by  Barratt  (1965). 
It  may  also  be  reexpressed  as  an  integral  with  respect  to  the 
wave  number  k  ,  through  the  connecting  relationship  given  by 
Eg .  (3.20).  The  result  is 


vC 

R  -  k  0 

k  tanh(kd)  j"p  2  +  p  2  +  02  +  0  2*1  ^  . 

2irpg 

. 

A  -  k  o  tanh  (kd)  /k  L  e  °  e  °  -1 

5< 

i 

(3.24) 

Here  kj  is  given  by  the  solution  of 


or 


kj  =  k  3  tanh(kjd)  if  k0d  >  1  , 
k,  =  O  if  k0d  <  1  . 


(3.25) 


It  may  be  noted  in  passing,  that  we  could  obtain  the 
steady  state  value  c..  the  total  power,  Eq.  (3.15),  in  a  sim¬ 
ilar  manner.  In  that  case,  after  integrating  Eq.  (3.15) 
with  respect  to  time,  we  would  get  the  following  expression, 
instead  of  Eq.  (3.16): 
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W.„  = 


c 

’2^2py 


W  dw 


du 


Y* 


P  2 
e 


P  2 
o 


+  Q  2  +  Q 


i 

0 


o*} 


<  - 


sin  (y 


+  wcj  t  sin(Y  -  wc) 
+  wc  y  -  wc 


_t 

■ 


As  in  the  case  for  R  ,  the  only  contribution  to  WT  , 
as  t  -*■<*>,  is  given  by  the  second  sine  term.  And  this 
occurs  when  y  =  wc  .  It  is  clear  then,  as  the  steady- 
state  condition  is  achieved,  that 


WT  =  CR 


=  WR  ' 


as  required. 


3.5  -  TWO-DIMENflONAL  WAVE  RESISTANCE: - 

It  is  of  interest  to  study  the  wave  resistance  of  a  two- 
dimensional  band  of  pressure  since  this  shows  up  clearly  the 
interference  effects  of  the  transverse  wave  system. 

A  possible  method  of  solving  this  problem  would  be  to 
set  up  a  two-dimensional  analog  of  Sec,(2)  -  for  the  potential, 
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,  -  and  then  solve  this.  However,  it  is  simpler  to  con¬ 
sider  a  pressure  distribution  independent  of  y  ,  and  then 
let.  the  beam  approach  infinity.  Thus  we  assume  that 


(x) 


cos 

sin 


(wx1* 


L 


(3.23) 


L  being  the  length  of  the  distribution. 

Let  us  now  consider  the  u  integral  in  Eq.  (3.13)  for 
the  wave  resistance.  We  have,  for  the  wave  resistance  per 


unit  width, 
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Tli  is  integral  is  given,  for  example,  in  Giadshteyn  and  Ryzhik 
(196'j,  p.  414,  3.741,  Formula  3).  Thus 


I  i 


71  COS 


/gwtanh (wd)  •  (t 


We  shall  now  examine  the  second  of  the  u  integrals, 
!;■  .  It  is  easily  seen  that 


I  1  < 

i2  I  _ 


bu‘ 


du  =  2//b 


l//b 


Thus,  as  b  -*•  ,  I2  -  0  . 

The  formula  for  the  two-dimensional  wave  resistance  may 
now  be  written  a s  the  following  double  integral: 


t 

r 


R 


2D  ’i  og 


CO 

f 


c  (  t  )  d  I 


dw  w 2  ♦ 


P2D2  +  Q2D2  IX 


■} 


;  *  < 


]-cos[. 


ycos  v gw • tanh (wd) • ( t  -  t)  -cos  w(s(t)  -  s(t))  .  (3.29) 
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3,6  -  TWO-DIMENSIONAL  TOTAL  POWER : - 

The  procedure  for  finding  the  2  uniting  expression  for 
the  total  power,  Eq.  (3.15),  in  the  case  of  two-dimensional 
flow,  is  identical  to  that  given  in  Sec.  (3.5).  The  result 
is 


W, 


T2D 


irpg 


c ( t )  di 


dw  w/gw*  tanh  (wd)  +  ^2D  ]' 


xsin  [/g  w*  tanh  (wd)  •  (t  -  x ) J* sin£w(s ( t )  -  s(x))J 


(3.30) 


3.7  -  TWO-DIMENSIONAL  STEADY-STATE  WAVE  RESISTANCE 

We  now  derive  the  two-dimensional  limit  of  Eq.  (3.22). 
The  procedure  followed  is  similar  to  that  in  Sec.  (3.5)  and 
the  result  is 
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If  k0d  ■  i 
there  is  no  real 
tancc  is  zero. 


(the  supercritical  speed  condition),  then 
solution  of  Eq.  (3.32),  and  the  wave  resis- 


4  -  STEADY-STATE  RESULTS 

4.1  -  PRESSURE  DISTRIBUTION  USED:- 

As  mentioned  in  Sec.  {1) ,  it  was  felt  that  it  would  be 
necessary  to  allow  the  cushion  pressure  to  drop  to  zero  wi¬ 
thin  a  finite  distance  at  the  periphery  of  the  craft.  The 
actual  details  of  this  fall-off  are  discussed  by  Hoaben  (1966b) 
and  Alexander  (1967)  for  a  periperal  jet  ACV,  and  by  Jones 
(1966)  for  a  plenum  chamber  machine. 

However,  it  is  unlikely  that  the  precise  way  the  pres¬ 
sure  drops  to  zero  is  crucial.  Rather,  the  essential  para¬ 
meter  is  a  measure  of  the  fall-off  distance.  The  restric¬ 
tion  of  planform  shape  to  a  rectangular  one  was  made,  since 
the  smoothing  effects  of  an  elliptical  design  were  studied 
by  Barratt  and  Newman.  The  general  pressure  distribution 
given  below  allows  the  fall-off  distances  to  be  individually 
varied: 

p(x,y)  =  po  £tanh  a(x  +  a)  -  tanh  a  (x  -  a)  Jx 

xftanh  8(y  +  b)  -  tanh  6  (y  -  b)l  .  (4.1) 


Here  po  is  the  nominal  pressure,  a  is  the  half-length 
and  b  is  the  half-beam.  The  smoothness  of  the  cushion 
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pressure*  fall-off  at  the  edges  is  a  function  of  the  fall-off 
parameters  t  and  .  A  large  value  of  a  implies  a 

rapid  drop  at  the  bow  and  stern  edges.  On  the  other  hand, 
if  3  is  large,  then  the  cushion  pressure  decays  rapidly  at 
the  sides.  This  function  is  shown  in  Fig.  2. 

As  a  particular  case,  we  may  consider  the  limit  as 
8  .  Then  the  pressure  is  given  by 

p  =  -^  p  o  £tanh  a(x  +  a)  -  tanh  a(x  -  a) 

for  -b  <  x  <  b  (4.2) 

=  0  otherwise. 

This  could  represent  a  so-called  sidewall  air-cushion  vehicle 
or  captured  air  bubble  (CAB) .  In  these  craft,  the  cushion 
is  physically  restrained  at  the  sides  so  that  the  pressure 
falls  abruptly. 

Another  particular  case  occurs  when  both  a  ♦  °°  and 
B  -*•  .  Then 


p=po  for  -a<x<a 

and  -b  <  y  <  b  (4.3) 

=  0  otherwise. 

This  case  of  a  constant  pressure  acting  over  a  rectangular 
area  is  one  commonly  used  by  previous  workers. 
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It  is  advantageous  to  obtain  analytic  expressions  for 
the  Kochin  functions,  Eq.  (3.14).  Due  to  symmetry  in  the 
chosen  expression  for  p  ,  it  is  clear  that 


Furthermore,  Eq.  (4.1)  allows  us  to  compute  P£  as  the  pro¬ 
duct  of  two  separate,  but  similar,  integrals  over  x  and  v, 
respectively.  The  contour  used  for  evaluating  the  integral 
is  shown  in  Fig.  3.  We  write  the  x  integral  as 


Ix  =  Real 


j^tanh 


a(x  +  a)  -  tanh  a(x  -  a)  e 


] 


1WX 


dx 


(4.5) 


Because  of  the  finite  lengths  of  the  paths  of  integration 
for  Ii  and  I3  in  Fig.  3,  it  is  clear  that  as  M  -+  ^  , 

1 1  “  1 3  =  0  .  Also,  along  C2  ,  points  corresponding  to  x 
on  C  are  given  by  x  +  in/a  .  It  can  then  be  shown  from 
Eq.  (4.5')  that 

I2  =  -  e_TTw/c<  I  .  (4.6) 

There  are  two  simple  poles  lying  inside  the  contour, 
and  after  applying  the  usual  residue  theorem,  and  some  sim¬ 
plification,  it  is  found  that 
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I 

x 


s  i  n  ( aw ) 

«  s  inh  ( nw/2  0 


(4.7) 


The 

tion,  so 


y  integral  in  P  may  be  written  down  by  inspec- 
that  the  final  result  for  the  Kochin  function  is 


_  r  sin  (aw)  tt  sin(bu) 

e  ^  0  a  sinh  (  nw/2  3  sinh(nu/2R) 


(4.8) 


has  well  defined  limits  for  the  cases  specified  by 
Eqs.  (4.2)  and  (4.3) . 

Finally,  by  integrating  the  pressure  distribution,  one 
obtains  the  weight  of  the  ACV : 

W  =  4  p  0  ab  .  (4.9) 


4.2  -  NONDIMENSIONAL  COEFFICIENTS 


The  wave  resistance  coefficient  used  is  defined  by 


R 

c 


R  aga 
w’  Po  ' 


and  the  power  coefficient  by 


W 

c 


W  c  >  ga 
W’  Pj 


(4.10) 


(4.11) 
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4.3  -  RESULTS 

The  programs  used  for  computing  the  steady-state  resis¬ 
tance  are  documented  by  Doctors  (1970,  Sec.  (10.1)). 

Some  curves  for  the  wave  resistance  are  shown  in  Figs.  4 
to  7  inclusive.  In  these  diagrams,  the  variable  used  for  the 
abscissa  is  A  =  1/(2F  )  .  This  results  in  the  low  speed 
interference  oscillations  having  a  period  approaching  a  con¬ 
stant  value  of  2*i  . 

With  regard  to  the  ordinate,  an  arcsinh  scale  is  oc¬ 
casionally  used  to  improve  the  presentation  of  the  results. 
This  transformation  is  defined  by 

yp  =  C i  arsinh(C2y)  ,  (4.12) 

where  y  is  the  number  to  be  plotted  and  yp  is  the  distance 
on  the  figure.  Ci  and  C2  are  constants. 

In  Fig.  4a,  the  eff?ct  of  smoothing  on  a  two-dimensional 
pressure  band,  that  is,  reducing  the  value  of  the  parameter 
aa  ,  is  illustrated.  For  aa  =  "  ,  we  have  a  sharp  pressure 
distribution  and  the  amplitude  of  the  oscillations  ^f  the 
wave  resistance  remains  unchanged  with  increasing  A  (de¬ 
creasing  speed) .  The  theory  predicts  an  infinite  number  of 
these  oscillations  between  zero  speed  and  any  finite  velocity. 
The  curve  is  a  sine  wave  and  this  result  was  obtained  by  Lamb 
(1932,  p.  403) . 

For  finite  values  of  the  pressure  fall-off  parameter, 
aa  ,  the  oscillations  die  out  with  decreasing  speed.  If 
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.a  =  5  ,  then  there  are  only  three  major  huisps.  The  humps 
and  hollows  are  caused  by  the  imerference  of  the  bow  and 
store  five  systems,  which  are  purely  transverse  in  this  ex - 
ample.  The  interference  is  not  affected  by  the  smoothing. 
However,  the  smoothing  reduces  the  amplitude  of  the  waves 
generated  ?y  the  bow  and  stem. 

Fig.  4b  shows  the  effect  of  smoothing  a  two-dimensional 

pressure  band  for  a  finite  depth.  The  resistance  is  zero 

for  supercritical  sneeds  (when  F,  >  )  )  because  free  waves 
*  *  a 

cannot  travel  faster  than  the  critical  speed. 

Three  cases  are  illustrated  ir.  rig.  5a.  Case  1,  with 
ja  =  a  -  *  ,  represents  a  sharp  pressure  distribution  with 
a  bean  to  length  ratio  of  0.5  traveling  over  deep  water. 

The  unrealistic  oscillations  obtained  by  previous  workers  are 
confirmed.  It  should  be  noted  that  for  low  speeds  the  amp¬ 
litude  of  the  oscillations  asymptotically  approaches  unity  - 
the  sane  as  for  the  tvo-dirensional  case.  Thus  the  trans¬ 
verse  waves  become  relatively  core  important  in  this  speed 
range.  Case  3,  with  ia  =  5a  -  3  ,  has  smoothing  applied 
on  u_ _  four  edges,  ar.d  once  again,  only  .about  three  cr  four 
humps  occur.  KSsiie  typical  values  of  aa  and  .-a  would 
have  tc  be  found  froo  measured  pressure  dis tributior.s ,  this 
curve  could  represent  an  actual  ACV.  Finally,  Case  2  only 
nas  smoothing  at  the  bow  and  stern  -  equivalent  to  a  side¬ 
wall  ACV.  The  result  is  almost  the  sane  as  for  Case  3, 
shoving  that  the  wave  pattern  is  essentially  produced  by  the 
fore  and  aft  portions  of  the  cushion  -  and  not  the  sides. 
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The  effect  of  smoothing  equally  on  all  four  edges  is  al¬ 
so  displayed  in  Fig.  5b.  The  sharply  euged  case  is  compared 
with  those  of  three  different  amounts  of  smoothing.  The  uia 
gram  shows  that  for  -_ia  =  Sa  >  5  an  unnaturally  large  num¬ 
ber  of  lev  speed  humps  and  hollows  occur. 

Fig.  5c  shows  the  same  three  cases  o:  r:- .  5a  L--t  for 
a  finite  depth.  The  chief  difference  now  is  that  th<--  main 
hump  is  pushed  tc  the  right  and  occurs  r.ear  the  critical  dept 
Froude  number.  The  low  speed  humps  are  only  .marginal  .y  af¬ 
fected. 

The  effect  of  varying  the  beam  to  length  ratio  is  de¬ 
picted  in  Fig.  6.  Fig.  6e  is  for  deco  water.  The  general 
effect  is  an  increase  in  the  maxima  of  the  resistance  coef¬ 
ficient  for  increasing  beam,  since  the  transverse  waves  as¬ 
sume  greater  importance  as  the  two-dimensional  case  is  ap¬ 
proached.  At  the  same  time,  the  minima  are  reduced  for  the 
same  reason.  A  secondary  effect  is  a  shift  in  the  location 
of  the  oscillations.  As  the  beam  increases,  this  shift  is 

to  the  right  (i.e,,  tc  lower  Froude  numbers).  In  the  two- 
dimensional  limit  the  hollows  occur  precisely  at  values  o: 

A  -  n-  ,  where  n  is  an  integer. 

The  result  of  varying  the  bear,  is  again  shewn  in  ~ig. 
but  now  for  a  finite  depth.  Similar  effects  are  shewn  in 
the  low  speed  range,  away  from  the  neighborhood  of  the  cri¬ 
tical  speed.  However,  the  position  of  the  critical  speed 
hump  is  hardly  affected.  For  finite  beams  it  occurs  at  a 
slightly  lower  speed  than  for  the  two-dimensional  case. 
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Fig.  7  illustrates  the  result  of  varying  tine  depth  of 
water.  Each  ot  tne  five  diagrams  is  for  a  different  beam 
tc  length  ratio.  They  show  how  the  critical  depth  hump 
resistance  increases  with  decreasing  depth.  For  the  two- 
dimensional  case,  this  peak  resistance  varies  as  the  inverse 
square  of  the  depth.  However,  for  finite  bean  to  length 
ratios,  this  maximum  drag  uoes  not  increase  nearly  as  rapidly 
in  shallow  water. 

Another  aspect  brought  out  by  these  diagrams  is  the  rate 
at  which  the  resistance  curves  for  different  depths,  but  the 
sane  bean,  converge  at  large  values  of  A  .  Thus,  for 
values  of  f.  greater  than  3  ,  the  case  of  d/a  =  1  is  in¬ 
distinguishable  frora  the  deep  water  case.  And  at  d/'a  =  0.5, 
the  resistance  curve  is  essentially  the  sane  as  that  for  any 
greater  depth,  provided  that  A  is  greater  than  6. 


5  -  CALCULATIONS  FOR  ACCULLRAYilD 


otio:; 


5.1  -  TWO-DIMKHSIONAL  RLSULYS : - 

Sample  results  for  accelerates  ~et:er.  are  s :vew:i  in 
Figs.  8  to  18.  All  cut  Fig.  9  were  computer  by  programs  lis 
tea  by  Doctors  (1970,  Sec.  (1C.2)).  Aitnouch  f'.e^e  pronr.vs 
can  handle  general  acceleration  patterns,  only  calculations 
corresponding  to  a  constant  acceleration  are  presented  sere. 

Various  two-dimensicnai  results  are  given  in  Figs.  8 
to  12.  Fig.  8  shows  the  effect  of  varying  the  level  of  ac¬ 
celeration  for  a  smooth  pressure  band  over  ceep  water. 

There  are  only  two  humps  displayed  when  the  acceleration  is 
nen-zero.  The  third  and  higher  order  humps  have  been  smooth 
ed  out  by  the  unsteady  motion.  Furthermore,  with  increased 
acceleration  up  to  0.2g  ,  the  second  last  hump  is  also  prac¬ 
tically  lost.  The  last  or  major  hump  is  also  somewhat  re¬ 
duced,  but  evidently  any  significant  reduction  can  only  be 
achieved  by  an  application  of  an  unnaturally  high  accelera¬ 
tion. 

The  location  of  the  humps  is  also  affected  by  the  ac¬ 
celeration.  With  increasing  acceleration,  the  oscillations 
are  delayed  to  higher  Froude  numbers. 

The  case  of  a  sharp  two-dimensional  pressure  band  is 
shown  in  Fig.  9.  A  complete  curve  could  not  be  computed  be¬ 
cause  of  the  (apparently)  infinitely  many  oscillations  that 
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occur  in  the  recjion  where  A  is  just  less  than  5  (for  the 
case  of  c''a  =  0.05).  The  last  two  humps,  and  the  last 
hollow,  correspond  in  position  with  those  displayed  for  a 
smooth  pressure  band  in  Fig.  3.  Thus,  the  smoothing  does 
not  affect  the  positions  of  the  oscillations  -  as  was  also 
evident,  for  steady  motion. 

The  resistance  coefficient  is  constant  and  equal  to 
unity  for  A  >  5  .  The  significance  of  this  point  lies  in 
the  fact  that  the  ACV  has  traveled  one  craft  length  when 
A  =  5  ,  at  an  acceleration  of  0.05g  . 

From  geometrical  considerations,  it  is  clear  that  there 
can  be  no  interference  effects  between  the  bew  and  stern 
transverse  waves  until  the  craft  has  moved  this  distance. 

Up  to  this  point  in  time,  the  bow  and  stern  individually  pro¬ 
duce  a  wave  train  of  varying  wavelength  (depending  on  the  in¬ 
stantaneous  velocity),  but  of  constant  amplitude.  bet  us 
call  this  amplitude  j;0  .  Then,  before  the  band  has  moved 
one  length,  there  are  two  wave  trains  produced,  each  con¬ 
taining  energy  per  unit  area  proportional  to  Co2  •  The 
wave  resistance  of  the  two  waves  in  this  speed  range  is  then 
proportional  to  2 r02  ,  and  is  constant. 

When  the  band  has  traveled  one  length,  the  stern  waves 
of  finite  wavelength  start  to  run  into  the  bow  system  that 
was  produced  at  the  start  of  the  motion.  These  latter  waves 
have  a  vanishingly  small  length.  As  a  result,  there  is  an 
infinite  accumulation  of  the  interference  oscillations  in 
the  resistance  curve  just  to  the  left  of  A  =  5  . 
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However,  later  on  in  the  motion,  the  two  wave  systems 
have  a  relatively  more  similar  wavelength,  or  frequency,  and 
then  have  a  capability  of  combining  tc  give  a  maximum  amp¬ 
litude  of  almost  2^c  ,  or  a  minimum  amplitude  near  zero. 
Hence,  the  wave  resistance  can  fluctuate  between  (~0)?  and 
at  the  higher  speeds.  So  the  i  eak  interference 
resistance  is  almost  twice  as  high  as  the  constant  value  oc¬ 
cur  ing  before  the  band  has  traveled  one  length. 

For  a  general  acceleration  level,  the  point  where  tine 
oscillations  start  is  given  by 

A j  =  g/4c  .  (5.1) 

Hence,  for  higher  accelerations.  A!  is  smaller  and  the  os¬ 
cillations  ara  displaced  further  to  the  left.  This  explains 
the  general  effect  of  acceleration  on  the  location  of  the 
humps  and  hollows  in  Fig.  8. 

Resistance  curves  for  a  smooth  distribution  appear  again 
in  Fig.  10,  but  the  water  now  has  a  finite  depth.  The  re¬ 
duction  in  magnitude  of  the  peaks  is  even  more  marked.  In¬ 
deed,  wnile  the  maximum  steady-state  resistance  increases 
without  limit  in  shallower  water,  the  peak  resistance  coef¬ 
ficient  for  accelerated  motion  is  more  nearly  fixed.  Thus 
the  relative  reduction  during  unsteady  notion  is  accentuated. 

The  location  of  the  unsteady  finite  depth  humps  is  af¬ 
fected  in  the  same  way  as  in  deep  water,  being  delayed  to  a 
higher  speed. 
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The  most  striking  phenomenon,  however,  is  that  for  all 
the  finite  depth  cases  studied,  the  wave  resistance  becomes 
negative  beyond  a  certain  velocity.  The  resistance  then 
asymptotically  approaches  zero  from  below. 

Other  interesting  features  of  the  two-dimensional  case 
are  presented  in  '•"’igs.  il  and  12.  These  show  the  total 
power,  resistance  power  and  sustention  power  -  which  were 
discussed  in  Secs. (3.1)  and  (3.3).  The  total  power,  being 
the  sum  of  the  other  two,  represents  the  total  rate  at  which 
work  is  being  done  on  the  water.  It  was  computed  in  order 
to  see  if  it  remained  positive  when  the  resistance  became 
negative.  Fig.  11  shows  these  three  powers  as  a  function  of 
time,  for  different  acceleration  levels  over  deep  water.  It 

is  seen  that  a  hump  in  the  resistance  power  or  total  power 
curves  generally  corresponds  to  a  hollow  in  the  sustention 
power  curve.  It  may  also  be  noted  that  the  sustention 
power  is  relatively  small  compared  with  the  two  other  quan¬ 
tities. 

Passing  now  to  the  case  of  finite  depth,  Fig.  12,  we 
notice  that  the  total  power  also  becomes  negative,  and  at 
approximately  the  same  point  in  time  as  the  resistance  chan¬ 
ges  iign.  In  shallower  water,  these  two  points  approach 
each  other.  At  a  depth  to  half-length  ratio  of  0.25,  the 
points  where  the  total  power  and  resistance  power  become 
negative  are  indistinguishable. 

In  interpreting  the  occurrence  of  negative  total  power, 
it  should  be  borne  in  mind  that  the  total  energy  input  to  the 
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water,  integrated  from  zero  time,  is  always  positive.  How¬ 
ever,  at  very  high  speeds  the  ACV  begins  to  recover  from  the 
wave  pattern  some  of  the  energy  previously  expended  in  form¬ 
ing  it.  The  location  of  the  negative  hump  in  the  resistance 
curve  is  further  discussed  in  Sec.  (6.6). 


5.2  -  THREE-DIMENSIONAL  RESULTS :- 

Some  computed  cases  for  a  three-dimensional  pressure 
distribution  are  shown  in  Figs.  13  to  18. 

The  resistance  over  deep  water  is  shown  in  Fig.  13. 

The  similarity  to  the  two-dimensional  analog  in  Fig.  8  is  to 
be  noted.  The  three-dimensional  resistance  is  generally 
about  half  that  of  the  two-dimensional  case  (for  this  beam 
to  length  ratio) .  Also  the  oscillations  occur  at  slightly 
higher  Froude  numbers,  due  to  the  additional  effect  of  the 
diverging  waves.  These  same  two  differences  display  them¬ 
selves  in  the  steady-state  results  in  Fig.  6. 

The  influence  of  acceleration  on  the  resistance  of  a 
three-dimensional  ACV  over  finite  depth  (Fig.  14)  is  seen  to 
be  remarkably  less  than  for  the  corresponding  two-dimensional 
situation.  Again,  it  must  be  remembered  that  the  steady- 
state  two-dimensional  results  showed  stronger  phenomena  than 
the  three-dimensional  ones.  Principally,  these  were  stronger 
interference  oscillations,  and  zero  resistance  beyond  the 
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critical  speed. 

Of  chief  interest  in  Fig.  14,  compared  vith  Fig.  10, 
is  the  lack  of  any  speed  range  in  which  the  resistance  is 
negative.  This  indicates  that  the  transverse  component  of 
the  wave  pattern  is  less  important  than  the  previous  steady- 
state  results  implied. 

The  three-dimensional  results,  nevertheless,  show  a 
substantial  reduction  in  the  peak  resistance  during  acceler¬ 
ated  motion,  particularly  for  the  shallower  situations. 

Beyond  the  last  hump,  the  curves  of  resistance  for  dif¬ 
ferent  levels  of  acceleration  converge  rapidly,  as  the  speed 
increases.  Beyond  A  =  0.5  (F  =  1.0)  the  acceleration  has 
little  effect. 

For  the  sake  of  completeness,  curves  of  total  power, 
resistance  power  and  sustention  power  are  given  in  Figs.  15 
and  16.  The  deep  water  case  (Fig.  15),  again,  is  similar 
to  the  two-dimensional  case  in  Fig.  11.  The  finite  depth 
results  in  Fig.  16  do  not  differ  considerably  from  the  deep 
water  results  of  the  previous  figure  -  in  contrast  to  the 
two-dimensional  results  discus&ed  before.  In  finite  depth 
the  three-dimensional  sustention  power  approaches  zero  some¬ 
what  faster  that  in  infinitely  deep  water. 

A  comparison  of  the  resistance  power  for  different  rates 
of  acceleration  is  made  in  Fig.  17  (deep  water)  and  Fig.  18 
(finite  depth).  For  accelerated  motion  over  deep  water, 
there  is  a  slight  increase  in  peak  power  to  overcome  resis¬ 
tance.  This  is  because  the  resistance  peaks  remain  essenti- 
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ally  the  same  in  magnitude,  but  are  shifted  to  a  hiyher  speed. 

At  a  depth  to  naif-  length  ratio  of  1.0,  the  resistance 
power  peaks  are  also  higher  in  accelerated  motion.  A  sig¬ 
nificant  reduction  in  the  maximum  resistance  power  occurs 
only  at  a  depth  to  half-length  ratio  of  0.25. 

Thus  we  have  a  situation  in  which  the  resistance  hump  is 
considerably  diminished  in  water  of  finite  depth.  However, 
the  new  location  of  the  hump  greatly  reduces  any  saving  in 
the  power  needed  to  overcome  it. 


*  -  ?CATirps  or  szxuxx  vvrr?. 


».  i  -  r^vc!£ : - 

It  is  wf!i  tr/vs  that  the  first-orJer  $ha!lcnr  vatet 
theory  sroisc*  i  '.cs*  realistic  results  for  ship  ^aw  resis¬ 
tance  than  the  I  i  wriied  theory  for  finite  death.  Thus 
’(irheli  ( 1 J1?  1  i»i  Tart  (I *f6i  shewed  that  it  predicted  zero 
resistance  ?or  a  to*y  travel*?^  at  less  thu  the  critical 
speed.  The  advantage  in  shali  ’jw  *4tc:  tbeopr  lies  in  the 
possibility  of  extending  it  to  cbtaia  hieber  order  approxi- 
■xat ions  to  the  flow,  veere  the  complexity  cf  the  finite  death 
results  vwIg  jreclsdt  this. 

The  theory  vas  deployed  is  this  case  to  predict  the 
l(y*it:a«j  of  the  negative  resistance  raorrp  displayed  ia  the  ms- 
steady  irfc-diweisivia  1  resales. 


6.2  ~  jgar/ATiog;  or  potsstials- 

The  procedure  is  essentially  the  sexe  as  for  the  case  cf 
finite  depth,  since  the  sane  set  of  equations,  naxelv  Ecs. 
(2.5),  [2. 7),  (?.9)  and  (2.115,  are  to  be  satisfied. 

srs  the  spirit  cf  the  shallow  water  anpronirration,  we 
shall  assure  ar.  asr^utotic  expansion  for  the  disturbance 
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S3 


•xstestial,  viz: 


L 

a=l 


\  9 


-  i) 


where  a1"  1  «■  o(s  >  -  it. 2) 

Purtheraore,  »»  shall  follow  Tsxrk  ard  siaoly  tar.e 

v  *  &({  I  .  ie.3) 

Sere  c  is  the  wraneter  de^oti-'w  the  shallowness  of  the 
vafer  sc  that  its  definition  cosld  be 

x,  y  =  0(1)  :  z  *  O(t)  ,  ifi.O 

x  „  y  ,  z  being  inside  the  fluid  recion. 

It  follows  froea  Eq.  (6.4J  that  dif ferent iation  can 
chance  the  order  of  saagnitude  of  a  quantity.  Thus 


Oil! 


(6.5) 


and 


0  ( e  " 1 )  - 


Ke 

Laplace 


nay  now  substitute  the  exoansicn, 
equation,  Eq.  *2.5).  Due  to  the 


F.q .  (6.1),  into 
assumptions  of 


the 
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Eqs.  (5.3)  and  i6.5),  we  say  separate  the  lowest  order  term 
in  Eq.  (2.5) .  Thus 

e(1)  *  0  .  (6.6) 

zz 

Siailarlv,  the  bed  condition,  Eq.  (2.11),  yields 
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[♦‘3,j  z-0  -  <(llt  +  c  -  0  '  (6-n) 

where  is  the  first  approximation  to  the  value  of  £  . 

The  last  boundary  condition  needed  is  the  dynamic  con¬ 
dition  on  the  free  surface,  Eq.  (2.9).  This  results  in 


[♦(1)fc  -  c  *(1)x]  Z=Q  +  p/p  +  g  C(1J  =  0  .  (6.12) 


The  last  three  equations  are  utilized  to  produce  the 

(l ) 

equation  to  be  solved  for  A  “  : 


-  2c  A 


-  c  A 


<X>  =  -  p 

X  p  x 


(6.13) 


The  method  of  solution  for  A'  '  is  similar  to  that 
used  in  Sec.  (2.2)  for  <j>  .  A  double  Fourier  transform 

(Eq.  (2.12))  with  respect  to  x  and  y  is  made,  and  this  is 
followed  by  the  substitution  given  by  Eq.  (2.19).  The  result 
of  these  two  operations  is 


Xtt  +  gd(w2  +  u2 )  x  =  — *p*exp(-iws  (t) )  . 


(6.14) 
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The*  similarity  between  the  previous  equation  and  the 
corresponding  one  for  finite  depth  should  be  observed.  The 
Laplace  transform  (Ea.  (2.13))  of  Eq.  (6.14)  is  taken,  and 
the  convolution  theorem  applied.  Then  the  double  inverse 
Fourier  transform  is  used  to  give  the  first  approximation  to 
the  shallow  water  potential: 


t  30  VC 


i'x,y,7.,t)  =  - - - 

f  f 

c'  dc*  :(-)  dT 

dw 

4  ti  2  o 

j 

■  l 

J 

j 

S’  0 


sin [/gd* k*  (t  -  t)]* 

vexp[i(w(x  -  x'  +  s  ( t)  -  s  ( t ) )  i  u(''  ~  y'))l,  '  1  :>) 

with  k2  =  w2  +  u2  . 
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6.3  -  HAVE  RESISTANCE 

The  first  aoproximation  to  the  resistance  is  just 


(1) 


(1) 


dS  . 


The  dynamic  condition  on  the  free  surface,  Eq.  (6.12),  is  now 
employed,  as  detailed  in  Sec.  (3.2).  The  final  result,  after 
simplification  is 


It  is  interesting  to  compare  this  result  with  that  for 
finite  depth,  namely  Eq.  (3.13).  The  shallow  water  result 
can  be  obtained  simply  by  letting  d  0  in  the  latter  to 
the  extent  that  tanh  (kd)  is  replaced  by  kd  . 

The  steady  state  resistance  can  be  obtained  either  as 
the  limit  of  Eq.  (6.16)  as  t  -+  »  (if  the  velocity  is  held 
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constant),  or  bv  allowing  d  to  approach  zero  in  Eq.  (3.24). 
The  result  follows: 


The  relation  between  the  transverse  and  longitudinal 
wave  numbers,  Eq.  (6.18),  shows  that  there  is  a  fixed  angle 
between  the  wavefronts  and  the  y  axis.  This  angle  is 
qiven  by 

0  =  arccos(l/Fd)  .  (6.1°) 

If  <  1  there  are  no  free  waves  satisfying  Fo.  (6.1^) 


and  the  wave  resistance  is  zero. 


6.4  -  ELEVATION  OF  THE  FREE  SURFACE 

The  free  surface  elevation  nay  be  found  from  substitut¬ 
ing  the  potential,  Eq.  (6.15),  into  the  dvnamic  condition, 

Eq.  (6.12).  Envployinq  the  standard  rules  for  differentiating 
an  integral  expression,  it  is  found  that 


c(1)  .  -  E_+  1 

s  pa  4  tt  2  og 


p'  dS' 


c(t)  dr 


dw  i  du 


"CO  “(X> 


w*cos [ /gd*k • (t  -  r ) ] x 

*sin[w(x  -  x'  +  s(t)  -  s  ( x ) )  +  u(y  -  y  * ) 3  .  (6.20) 


On  the  other  hand,  if  we  restrict  our  attention  to 
pressure  distributions  which  are  symmetric  about  the  x  and 
y  axes,  such  that 


P 

o 


=  0  , 


then  Eq.  (6.20)  may  be  simplified  to 
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6.5  -  STEADY-STATE  RESULTS 

A  comparison  of  shallow  water  and  finite  depth  theorv 
is  made  in  Fig.  19.  The  case  of  a  smooth  three-dimensional 
pressure  distribution  is  depicted  in  Fig.  19a.  It  is  seen 
that  a  depth  to  half-length  ratio  of  0.25  is  essentially 
shallow  if  1/F^2  *  0.5  . 

The  wave  resistance  of  a  sharp  three-dimensional  dis¬ 
tribution  is  shown  in  Fig.  19b.  The  convergence  of  the 
finite  depth  resistance  to  the  shallow  water  one  is  per- 
hans  not  so  ranid  her-'.  An  inter estinci  feature  of  the 
shallow  w  Ler  i esul 1  Is  the  knuckle  in  the  resistance  curve 

fi'ho  h  her  sudden  chances  in  slope  in  the 
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curve  are  due  to  the  fact  that  a  finite  number  of  points  have 
been  computed  to  define  it.)  For  this  pressure  distribution, 
it  is  possible  to  perforin  the  integration  in  Eq.  (6.17)  analyt¬ 
ically.  Using  Gradshteyn  and  Ryzhik  (1965,  p.  451,  3,828, 
Formula  9) ,  the  result  may  be  written  down  as 

F  2 

r  =  a-  ,  d 
c  b  2{  Fd2  -  1)  3/2 

V 

2(Fd2  -  1) 


where  F,  =  /a2/b2  +  1  . 

d  i 

Thus,  the  resistance  is  continuous  at  a  depth  Froude 
number  equal  to  Fdj  but  the  slope  is  not.  It  is  easilv 
verified  that  when  Fd  =  Fd  ,  the  wave  pattern  appears  as  in 
Fig.  20.  Apparently  there  is  a  reinforcement  between  the 
bow  and  stern  wave  systems  at  this  particular  speed,  causing 
the  observed  knuckle. 

In  Fig.  19a  a  slight  hump  is  also  discernible  at  this 
depth  Froude  number,  but  due  to  the  finite  values  of  the 
pressure  fall-off  parameters  the  slope  of  the  resistance 


if  Fd  >  Fd, 


(6.22) 


if  Fd  <  ?d,  ' 


curve  is  continuous. 
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6.6  -  AC CE LE RAT E D  MOT I ON ; - 

In  order  to  examine  the  phenomena  exhibited  in  the  two- 
dimensional  wave  resistance  curves  of  Sec.  (5.1),  we  shall 
consider  the  free  surface  elevation  due  to  the  motion  of  a 
pressure  band  qivon  by  Eq.  (4.2).  The  Kochin  function,  , 

(liven  by  Eq.  (4.8)  is  then  substituted  into  Ea.  (6.21). 

The  two-dimensional  limit  as  b  -*■  «>  is  obtained,  by  use  of 
Eq.  (3.18)  to  evaluate  the  u  integral  and  the  result  is 


'O  t 

4  |  dw  c  (  t  )  *w  •  cos  [ ,  qd  *w  •  ( t 

oq  pact  | 

*  I  | 

0  0 


T)  )  v 


'sinlwlx  +  s(t)  -  sh)11'sinMi!w/L)  dl 


(6.23) 


We  now  apply  the  formulas  for  the  products  of  sines  and 
cosines,  twice,  to  give 


r.(1)  =  -  E-  +  l  I  i 

og  ifj-h 


<«  t 


i  •  P-°-  |  dw  |  c(t)  >wx 

4ogct  1  ' 


0  0 


cos  [w  ( 1/gd  •  (t  r) _ 4-  j  (x  4-  s(t)  -  s(t)  -  ia)  )  ] 

sinh  (irw/2a) 


(6.24; 


sinus  <_  f  the  impropriate 


total  nf  our  eims  in  the  double 
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summation. 

The  w  integral  may  now  be  done  usinq  Gradshteyn  and 
Ryzhik  (1965,  p.  511,  4 . Ill ,  Formula  2): 


l  I  i 

i, j=±l 


t 

■ 

0 


c  (t)  x 


xsech2  {a  ( /gd*  (t  -  x)  +  j  (x  s(t)  -  s(t)  -  ia)  )  J  dx  .  (6.25) 


Further  restricting  our  interest  to  a  sham  pressure 
band,  so  that  a  -*  °°  ,  it  is  seen  that  the  integrand  is  zero 
except  when 


v^gd*  (t  -  t)  +  j  (x  +  s(t)  -  s(x)  -  ia)  =  0  .  (6.26) 


Let  us  say  that  this  occurs  when  t  =  t,  and  the  velocit"  of 
the  band  is  c(ti)  =  Cj  ,  while  the  distance  traveled  is 
s(tj)  =  sj  .  In  the  neighborhood  of  x  =  tj  ,  we  put 


T  =  t  j  -i-  T  '  , 


where  t*  is  small  as  a  -*•  «»  .  The  arnument  of  the  hvoer- 
bolic  secant  in  Eg.  (6.25)  may  then  be  reexpressed  in  the 
region  where  t  =  t  .  It  is  just 

t"  =  -  a  (/gd  +  jcj )  t  '  . 


(6.25)  to  nivc 


f.  4 


i-'il t  i  tut  ion  is  now  mado  in  i’n , 


i  ,  !  •  •  i  4  pa  (  I  gTT  1  jc  i  ) 


"•  sgn  ( ''yd  jc  ,  ) 

j  soch-  r"  dx"  .  (6.27) 

( 

'*>  ’  s an  ( »'gd-t  j  c  j  ) 


Tn  this  integral ,  the*  unoor  limit  should  bo  replaced  bv  zero 
*  !'  t  I  ■■  t  ,  and  similarly  for  the  1  -wer  Limit.,  if  t  -  =  0  . 
Wo  then  obtain 


(I)  =  pa 

pq 


Hero 


an  : 


-  H  ( a  -  !  x  |  ) 


i »  j  =  ± 1  I  1  +  j -F 


1  if  0  <:  t,  <  t 


4  if 


0  or  t  i 


0  if  h  i  <  0  or  t  ]  '•  t  , 


'VISU1:'  S' 


i me*"  ■  on  . 


((■  .23) 


( r, .  2  f  i ) 


AVAILABLE  COPY 


6S 


The  wave  resistance  of  this  sharn  pressure  band  *av  b*» 
derived  from  the  forsula: 


0(1) 

?  2D  "  ®* 


i=-i  *-  J  x  -* 


The  resistance  coefficient  then  follcvs: 


The  shallow  water  resistance  is  clotted  as  s  function 
of  tine  in  Fio.  21.  This  is  for  a  constant  acceleration . 

At  one  point  the  resistance  becomes  infinite,  and  at  another 
it  becomes  negatively  infinite.  Ti.is  second  ceak  corres- 
oonds  closely  with  the  negative  humn  Dredicted  by  the  finite 
depth  theory  for  the  same  deoth  and  acceleration  (Fin.  12e)  . 
Apparently  the  positive  hump  does  not  show  uo  in  finite  deoth 
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theorv.  It  should  ho  noted  that  the  positive  humo  in  the 
shallow  water  *■  heorv  is  snoother  than  the  neoative  one  in 
the  sense  that  the  latter  has  an  infinite  discontinuity  in 
function  value. 

Vi  i.  22  shows  *  he  shane  of  the  water  surface  at  dif¬ 
ferent  sta«ies  of  the  motion  of  Pin.  21.  Initially,  the  de¬ 
pression  is  due  only  to  the  static  pressure  under  the  pres¬ 
sure  band.  As  time  prooresses  a  pile  of  water  lids  up 
ahead  of  the  bow  and  this  develops  into  an  infinite  peak  at 
the  critical  depth  Froude  number  (r  =  1).  At  the  same  tine, 
a  negative  peak  develops  iust  inside  the  The  p-iir  of 

peaks  causes  the  positiv  •  hump  sho’-ai  in  ;  > .  21. 

Bevond  the  critical  speed  these  peaks  cannot  keen  up 
with  the  ACV ,  and  are  seen  to  move  back  relative  to  it.  At 
a  second  critical  point  m  time,  the  bow  peak  passes  un  ler 
the  stern  of  the  pressure  band  and  the  resistance  sud  ienlv 
drops  to  minus  infinity.  Thereafter,  the  water  surface  in 
the  vicinity  of  the  ACV  levels  out,  and  the  wave  resistance 
approaches  zero. 

An  inspection  of  Eq .  (6.30)  shows  that  the  peaks  in  the 
resistance  curve  occur  when  i  -  -1  and  F  ,  =  1  .  Look- 

di 

inn  now  at  Eq.  (6.31),  this  can  happen  when: 


Tiiis  ■  li  :  L  ..  t;  other  ~ .  •  s.ibilitv 

is  • 
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ft  -  -i :  v/gd*  (t  -  ti)  -  (s  -  s  i  4  2a)  -  0  . 

This  case  corresponds  to  the  negative  peak.  Of  the  two  con¬ 
ditions  represented  by  the  latter  equation,  the  only  one  that 
satisfies  the  relationship  0  <  tj  <  t  is 


Eq.  (6.32)  gives  the  value  of  t  when  the  soike  of 
water  generated  at  the  bow  passes  under  the  stern  of  the 
pressure  band. 

If  we  now  return  to  the  special  case  of  a  constant  ac¬ 
celeration  c  ,  Eq.  (6.32)  becomes  a  quadratic  in  t  .  The 
depth  Froude  number  which  locates  the  negative  humD  is  then 
given  by  the  formula: 


(6.33) 
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Kq.  (6.33)  is  plotted  in  Fiq.  23,  toqether  with  some 
points  represontinq  the  location  of  the  necrative  hump  for  a 
finite  depth.  As  expected,  the  agreement  is  better  for 
shallower  depths  and  sharper  pressure  distributions. 


7  -  THE  INVERSE  PROBLEM 


7.1  -  PROBLEM  STATEMENT 

An  object  of  interest  in  the  design  of  an  ACV  is  to 
know  the  attainable  acceleration  pattern  as  a  function  of  the 
characteristics  of  the  propulsor.  Hence,  other  items  of 
concern,  such  as  the  acceleration  margin  at  the  harm  speed, 
may  be  deduced. 

In  order  to  make  the  study  realistic,  one  must  include 
all  the  components  of  resistance.  These  were  described  in 
Sec.  (1.1).  We  use  Newton's  law  to  obtain 

T  -  D.  -  D..  -  R  =  me  ,  (7.1) 

A  M 

where  T  is  the  thrust  available,  is  the  aerodvnamic  or 
profile  drag,  DM  is  the  momentum  dracr  and  m  is  the  mass 
of  the  craft.  The  problem  definition  could  be  made  more 
elaborate  by  including,  for  example,  an  estimate  of  the  water 
contact  resistance.  It  is  emphasized  here  that  this  break¬ 
down  into  drag  components  is  only  an  idealization  of  the  true 
situation  in  which  there  are  interactions  among  them. 

In  Eq.  (7.1),  the  terms  on  the  left-hand  side  are  func¬ 
tions  of  speed,  or  acceleration  pattern  -  which  are  unknown. 
We  shall  assume  the  following  forms  of  these  resistance  com¬ 
ponents  : 
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n,  - 


2  CD  AF  r  A 


(7.2) 


and 


nM  = 

M 


mA  C  ' 


(7.3) 


where  is  the  aerodvnamic  drag  coefficient  of  the  craft 

based  on  A„  ,  the  frontal  area,  v'hile  cA  is  the  air  den- 
sity  and  is  the  mass  flow  rate  into  the  air  cushion. 

The  program  used  to  solve  Eq.  (7.1)  is  described  by  P  >clors 
(1970,  Sec.  (10.3)). 

Let  us  now  limit  the  numerical  investigation  to  a  single 
craft.  The  data  are  aiven  meaninq  bv  relating  them  to  a 
particular  (dimensional)  machine ,  below.  The  results  will, 
also  apply  to  any  scale  model  of  this  craft. 


DIMENSIONAL  VARIABLE 

DIMENSIONLESS 

’LAriAPLE 

P 

= 

1.94 

slugs/ft . 3 

(base  unit) 

q 

= 

32.2 

ft . /sec . 2 

(base  unit) 

a 

- 

20 

ft. 

(base  unit) 

b 

= 

10 

ft. 

b/a  = 

o  .  5 

a 

= 

0.25 

ft.-1 

a  ci  ~ 

5 

P 

= 

0.25 

ft.-1 

3a 

1 

5  1 

w 

- 

25,000 

lbf . 

pga/po 

/n 

o 

\ 

- 

0.002425 

slugs. /ft. 3 

■  a  1 

"  -  ■  2 '  - 

r.  -  ’60  ■  .  1  v  *.m  ■  '  )(j'; 

;  j  1 1  * 


".  i  34 


i  Lugs/ sec . 


0.0075 
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7.2  -  RESULTS  FOR  UNSTEADY  MOTION 

Fig.  24  shows  the  effect  of  applying  three  different 
constant  thrust  levels  to  the  ACV.  At  the  hump  soeed  the 

wave  resistance  is  of  the  sane  order  as  the  two  other  draq 
components  combined.  Only  one  other  minor  hump  is  dis¬ 
played,  as  was  also  the  case  for  a  constant  acceleration 
level.  Thus  the  character  of  the  wave  resistance  curve  is 
not  grossly  affected  by  the  nature  of  the  acceleration  pat¬ 
tern  -  provided  it  is  reasonably  smooth.  At  the  lowest 
thrust  level  (0.06  of  the  craft  weight)  the  acceleration  mar¬ 
gin  at  the  hump  is  onlv  about  0.02q  which  could  be  inade¬ 
quate  for  a  practical  machine. 

The  same  craft  running  over  finite  depth  is  presented 
in  Fig.  25.  Similar  qeneral  trends  are  displayed.  However, 
the  thrust  margin  at  the  hump  speed  is  slightlv  reduced,  but 
thereafter  the  wave  resistance  drops  somewhat  faster  than  in 
deep  water, 

The  corresponding  velocity  patterns  are  shown  in  Fias. 

26  and  27.  The  curves  are  quite  smooth  despite  the  humps 
in  the  wave  resistance.  The  effect  of  depth  is  seen  to  be 
small  for  this  case  -  the  dron  in  velocity  in  the  reoion  of 
the  hump  due  to  finite  depth  is  less  than  10?,. 

Some  calculations  were  also  performed  usina  an  enaire- 
driven  airscrew  to  push  the  craft  rather  than  a  constant- 
thrust  device.  The  propeller  used  was  a  four-bladed  Clark  Y 
Section  screw  (Propeller  Reference  Number  5868-9) .  The 
characteristics  of  this  propeller  were  measured  by  Hartman 
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and  liinrntann  (1977).  The  diameter,  D  ,  was  10  ft.  and  the 
blades  had  a  nominal  angle  of  20°  at  the  three-quarter  indial 
ooint.  The  table  below  lists  the  thrust  and  toroue  eoeftic  i- 
ents,  and  K  ,  as  a  function  of  the  advance  ratio,  J  , 

where 

kt  =  t/*>an2dU  '  (7-4) 

K0  =  VdaN-d5  (7.5) 

and  J  =  c/MD  .  f  7 . 0 ) 

Here  Q  is  the  propeller  torque  and  N  its  speed  <n  >  e/o 
lutions  per  unit  time. 
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Fig.  28  shows  the  variation  of  drag  components  with  time, 
assuming  that  the  propeller  absorbs  a  power,  P  ,  given  by 

P  =  3. 804x10s  ft.lbf./sec.  (692  H.P.)  , 


or  P/ma  */2g  3//2  =  0.6  . 

The  results  are  very  similar  to  those  in  ^ic.  24  for  a  cons¬ 
tant  thrust.  This  is  because  the  propeller  advance  ratio  is 
still  small  at  the  hump  speed  -  and  up  to  this  speed 
and  Kq  are  almost  constant 

The  case  of  a  constant  propeller  speed  of  revolution  is 
treated  in  Fig.  29.  Fig.  29a  hows  the  resistance  components 
when  the  propeller  runs  at  a  speed  given  by 

N  =  25,38  r.p.s.  , 

or  N/a/g  =  20  . 

Again,  the  thrust  drops  only  slightly  in  the  speed  ranee  in¬ 
volved,  so  the  results  are  similar  to  those  for  a  constant 
thrust,  or  a  const'nt  encine  power.  In  Fies.  29b  and  29c 
the  craft  is  overloaded  by  50%  and  100%,  resoectivelv .  Onlv 
the  weight  (and  cushion  pressure)  were  increased,  and  all  the 
other  dimensional  data  were  held  fixed.  Since  the  wave 
drag  is  proportional  to  the  square  of  the  cushion  pressure, 
this  component  becomes  relatively  more  important  in  cases  of 
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overload.  For  the  100%  overload  situation,  the  machine  hare- 
lv  surpasses  the  hump. 

The  case  of  the  same  propeller  SDeed  of  revolution,  but 
over  finite  depth,  is  handled  in  Fig.  30.  For  a  normal  load 
(Fiq.  30a),  or  a  50%  overload  (Fig.  30b),  the  drag  components 
follow  quite  closely  the  corresponding  curves  for  deep  water 
given  in  Figs.  29a  and  29b,  respectively.  However,  for  a 
100%  overload  (Fiq.  30c)  the  situation  is  quite  different  in 
finite  depth. 

Here  the  craft  apparently  does  not  overcome  the  resis¬ 
tance  hump  in  the  time  shown.  At  first,  .it  seems  fo  onlv 
just  cross  the  (unsteady)  hump.  Then  the  wave  resistance 
rises  -  probably  because  the  almost  zero  acceleration  at  this 
time  allows  the  shallow  water  wave  pattern  to  build  up,  thus 
approaching  a  nearly  steady  state.  The  total  draq  (beina 
essentially  the  same  as  the  wave  resistance  now)  is  creator 

than  the  thrust  so  that  the  craft  decelerates  to  belo1-  the 

hump  speed.  Without  extending  the  calculation  to  a  deal¬ 
er  time  limit,  it  is  difficult  to  predict  whether  the  v  V 
would  cross  the  hump,  approach  a  steadv-state  sub-hump  speed, 
or  settle  dovm  into  a  repetitive  cvcle. 

The  velocity  pattern  for  the  above-mentioned  **hi  '  load¬ 
ings  is  shown  in  Mias .  31  and  32.  In  deep  water  (Fii.  11) 
it  is  seen  that  'while  a  50"  overload  ■  tolerable  for  'bis 
craft,  a  KKM  ove  •  low  7  •.•iv-.ld  or  ever  1  b.une  _.d  be  ,  nu  achiev¬ 
ed  in  o'.1 'lice.  i  ri  ■  ’  :  ' (Fie.  '  V  the  rr-  i  i  <n  of  de¬ 
coloration  f-v  IPO'  o  -|  i  :  mow  "u'i  e  c  1  u<  rl\  roi  o.i !  wd  . 


The  main  assumption  here  is  that  the  propeller  and  en¬ 
gine  characteristics  at  any  instant  of  an  unsteady  notion  can 
be  estimated  on  a  quasi-steady  basis  using  only  the  instan¬ 
taneous  craft  speed,  thus  ignoring  inertia  and  rime  history 
effects.  (In  fact,  the  only  quantity  treated  in  a  truly  un¬ 
steady  fashion  here  is  the  wave  resistance.)  The  calculated 
results,  retrospectively,  seen  to  confir.r.  the  validity  of  t..i 
assumption.  The  required  changes  in  prop? Her  speed  of  ro¬ 
tation  (under  the  constraint  of  constant  power)  and  t.nrust 
(under  the  constraint  of  constant  propeller  speed)  are  indeed 
quite  small.  Thus  the  system  of  the  propeller  and  engine  is 
effectively  operating  at  a  nearly  steady  state. 


7.3  -  QUASI-STEADY  RESULTS 

One  might  now  consider  the  effect  of  a  further  simpli¬ 
fication,  namely  the  assumption  that  the  wave  resistance  can 
also  be  treated  on  a  quasi-steady  basis.  In  other  words, 
the  wave  resistance  (as  well  as  the  other  drag  components) 
at  any  instant  of  the  unsteady  motion  is  assumed  to  be  in¬ 
dependent  of  acceleration  and  time  history  -  and  just  equal 
to  the  steady-state  value  at  the  instantaneous  velocity. 

This  technique  could  make  use  of  the  steady-state  results 
previously  calculated  in  Sec.  (4.3),  which  incidentally 
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require  saucn  less  vomputat lonu  1  effort. 


tno  dr^g  components  over  deep  water  for  tnree  different 
thrust  levels.  'i is  generally  little  difference  in  tie 
resultant  motion  oecause,  as  was  seen  in  S<  e.  (5.2),  t  ie  mag¬ 
nitude  of  the  hump  resistance  i:i  deep  water  is  only  slightly 
affected  by  this  range  of  acceleration.  Moreover,  toe  low- 
speed  unsteady  resistance  curve  tended  to  pass  through  the 
average  value  of  tae  oscillating  steady-state  resistance 
curve,  and  tne  effect  of  small  local  oscillations  ir.  asis- 
tanee  tends  to  be  averaged  out  in  the  resultant  velocity  pat¬ 
tern. 


The  effect  of  finite  depth  is  shown  in  Fig.  34.  The 
difference  compared  witn  the  truly  unsteady  calculation 
(Fig.  25)  is  more  discernible  row  -  because  of  the  higher 
steady-state  hump  encountered  in  the  quasi-steady  calcula¬ 
tions.  This  aspect  is  .more  clearly  seen  in  Figs.  i.  an:  ’6 
which  display  the  velocity  patterns  for  deep  water  end  finite 
depth,  -espectively .  For  comparison,  the  unsteady  curves 
are  also  shown.  In  deep  water  there  is  little  difference, 
but  for  d/a  =  0.5  ,  the  patterns  separate  slightly  -r- 
the  hump  speed  for  the  reasons  cited  cbc-ve. 

Toe  result-  of  employing  the  propeller  previous!"  des¬ 
cribed  (with  a  constant  speed  of  revolution)  is  shown  in 
r‘:  : .  i~  da-' ~  .  v  »:  i  :.  ;J  (  5)  .  A ;  was  n  >t  vf 

in  Sec .  .  j.)  ,  •  v  .*  Tv.‘l  1  m  a ••-‘u  v  •  u  c. 

to  a  cons*”  n.t  .  •  : .  *.  .  Bel  ••  •  n?  i  inp  sp-  1 ,  flic 
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propeller  produces  an  almost  constant  thrust.  In  deep  water 
the  quasi-steady  and  unsteady  results  (Fig.  29)  are  again 
similar. 

For  d/a  =  0.25  ,  the  motion  of  a  normally  loaded  craft 
compares  with  that  for  the  unsteady  calculation  (Fig.  30a). 
However,  for  a  50%  overload,  the  quasi* steady  calculation 
predicts  that  the  ACV  can  only  achieve  a  sub-hump  speed  while 
the  unsteady  theory  (Fig.  30o)  shows  tnat  the  ACV,  in  fact, 
has  a  practical  thrust  margin  to  surpass  the  hump  resistance. 
In  the  case  of  100%  overload  the  unsteady  result  (Fig.  30c) 
shows  that  the  machine  is  in  a  marginal  position  regarding  its 
ability  to  achieve  a  practical  cruising  speed.  The  quasi - 
steady  calculation,  or.  the  other  hand,  would  indicate  that 
the  craft  considerably  lacks  the  ability  to  cross  the  hump  - 
it  would  seen  to  have  only  about  half  of  the  thrust  required. 

The  velocity  patterns  for  a  constant  propeller  speed  of 
revolution  are  shown  in  Fig.  39  (deep  water)  and  Fig.  40 
(d/a  -  0,25).  These  clearly  indicate  tiie  narked  difference 
between  the  predictions  of  the  cuasi-steadv  and  unsteady  cal¬ 
culations  -  particularly  for  finite  depth. 


8  -  CONCLUDING  REMARKS 


8.1  -  CONCLUSIONS : - 

A  large  number  of  results  have  been  presented,  and  the 
reader  is  referred  to  the  appropriate  section  for  the  de¬ 
tails.  However,  some  general  remarks  should  now  be  made, 
and  the  first  of  these  concerns  the  degree  of  smoothing  re¬ 
quired  to  eliminate  the  low  speed  oscillations  in  the  steady- 
state  wave  resistance  curves. 

Deference  is  now  made  to  Fig.  2  showir.a  the  pressure 
distribution.  We  shi  11  take  2a  =  <u  since  this  parameter 
hardly  affects  the  resistance.  The  slope  of  the  free  sur¬ 
face  at  the  center  of  the  bow  (x  =  a,  y  =  0) ,  at  zero  speed, 
is  then  given  by 


x 


sech  2  (2na) 


For  r,a  ^  2  the  second  term  may  be  dropped  to  give 


r  =  SEl  . 
x  *  2ng 


(8.1) 


S'->pir  cnluul  at  i  u  i  v  toner  (  I  ;'<•{?)  .  “or  cxei  ••  i.e,  indic¬ 
ate  tii>'  t !  ■  -  s  .  -  •>>.'  i  1  !'■  t»> . ;  «*i  unit-.  tr  l  l  .  Tut'  ‘t'o , 

for  t.vpj  \!  1  C's,  the  rifio  r'ia/pn 


in  about  40. 
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Eg.  (8.1)  would  yield  aa  =80  as  a  general  value.  'T’ho  con¬ 
clusion  in  this  regard  then,  is  that  a  reasonable  amount  of 
smoothing  is  not  sufficient  by  itself  to  eliminate  the  low 
speed  humps  and  hollows.  Other  effects,  such  as  nonlinenritv 
and  real  fluid  properties,  may  have  to  be  invoiced,  to  explain 
the  discrepancy. 

It  is  evident,  however,  that  the  use  of  low  values  of 
aa  and  $a  (about  5)  does  improve  the  agreement  between  cal¬ 
culated  and  observed  steady-state  wave  resistance  curves,  es¬ 
pecially  at  lower  speeds.  Thus  the  use  of  such  values,  oven 
if  unnecessary  from  the  point  of  view  of  accurate  pressure 
modelling,  could  still  be  justified  if  aa  and  Pa  are  re¬ 
garded  as  semi -empirical  factors  introduced  to  effectively 
eliminate  the  unrealistic  linear  effects  at  low  Froudo  num¬ 
bers.  Tvo  advantages  wall  accrue. 

First,  the  inverse  problem  can  be  treated  with  more  con¬ 
fidence  since  it  requires  a  continuous  calculation  of  wave 
resistance  for  the  entire  speed  range  beginning  from  zero  and 
one  can  avoid  the  unrealistic  phenomena  which  would  occur  for 
a  sharp  pressure  distribution. 

Second,  the  computational  effort  involved  in  evaluating 
the  various  integrals  is  considerably  reduced  duo  to  the  ex¬ 
ponential  decay  of  the  integrands  at  infinitv  if  aa  and  -a 
are  finite  and  small.  In  fact  the  correspondinn  calculations 
for  a  sharp  pressure  distribution  would  be  formidable,  as  dis¬ 
cussed  in  Sec.  (8.2).  Thus  there  is  also  a  practical  and 
economical  advantage  in  the  use  of  these  adjustable  constants. 
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Most  of  the  figures  show  a  marked  similarity  between  the 
two-  and  three-dimensional  results.  Thus  the  transverse  wave 
system  plays  a  large  role  in  contributing  to  the  energy  radi¬ 
ation  for  the  usual  beam  to  length  ratios  of  ACVs .  On  the 
other  hand,  the  diverging  waves  are  sufficiently  important  to 
prevent  the  occurrence  of  negative  wave  resistance  durinn  ac¬ 
celerated  motion  of  the  craft  from  rest  in  finite  depth. 

The  effect  of  acceleration  level  on  the  wave  resistance 
is  quite  marked.  Normal  accelerations  of  up  to  0.1c  are 
sufficient  to  eliminate  most  of  the  low  speed  humps.  Also, 
the  high  resistance  peak  at  the  critical  dnpth  Froude  number 
is  limited  to  a  much  smaller  value  in  accelerated  motion. 

However  the  resistance  peak  is  displaced  to  a  higher 
speed  so  that  the  power  needed  to  overcome  it  is  hardly 
diminished,  except  in  very  shallow  water. 

On  the  other  hand,  the  practical  design  of  an  ACV  results 
in  the  propulsive  efficiencv  being  a  maximum  near  the  design 
speed.  At  low  speeds,  up  to  the  hump,  the  thrust  is  essenti¬ 
ally  constant.  Hence  the  critical  tiling  is  the  peak  wave 
resistance  -  and  not  the  peak  power  to  overcome  it.  This  was 
borne  out  by  a  comparison  of  the  unsteady  and  quasi-stendv 
theories  for  the  solution  of  the  inverse  problem  in  finite 
depth.  Under  certain  circumstances  the  ACV  could  easi.1” 
cross  the  hump,  even  though  the  quasi-steadv  theory  predicted 


o the r  ■  i so . 
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8.2  -  FUTURE  WORK;- 

It  would  be  worthwhile  to  run  some  experiments  to  veri¬ 
fy  the  calculations  for  accelerated  motion.  Either  the 
direct  or  inverse  case  could  be  tested  since  the  theory  is 
equally  valid  for  both.  Correction  for  the  additional  drag 
components  is  one  of  the  problems  to  be  contended  with. 

A  particular  case  that  should  be  examined  is  the  two- 
dimensional  one.  It  should  be  interesting  to  see  how  ac¬ 
curately  the  theory  predicts  the  phenomenon  of  negative  wave 
resistance . 

A  few  comments  should  also  be  addressed  to  the  problem 
of  computing  the  integral  in  Eq.  (3.13).  The  convergence  of 
this  integral  deteriorates  for  large  values  of  t  ,  a  and 
8  .  This  is  because  of  the  exponential  decays  in  the  Kochin 
function,  Eq.  (4.8).  The  Gauss-Laguerre  quadrature  rules 
used  are  based  on  these  decays.  Thus  the  majority  of  the 
resistance  curves  for  accelerated  motion,  given  in  the  figures, 
employed  a  rule  of  order  2048  for  the  w  integral,  and  a  rule 
of  order  16  for  the  u  integral.  However,  an  estimation  of 
the  truncation  error  showed  that,  in  fact,  only  the  first  102 
and  9  points,  respectively,  were  needed  for  a  three-figure 
accuracy.  Generally,  a  129-point  trapezoidal  rule  was  used 
in  the  t  integral. 

This  combination  of  rules  was  found  to  be  the  best,  al¬ 
though  it  would  be  useful  to  investigate  other  numerical  quad¬ 
rature  schemes. 


Because  of  the  limitations  of  the  numerical  method 
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specially  for  large  t  -  an  alternative  technique  of  evalu¬ 
ating  the  integral  should  be  considered.  For  values  of  the 
velocity  somewhat  greater  than  the  main  hump  speed,  it  should 
be  possible  to  construct  an  asymptotic  expression  for  the 
wave  resistance  (for  large  t) .  For  the  inverse  problem, 
the  asymptotic  formula  would  have  to  be  based  on  an  extra¬ 
polated  velocity  pattern  (which  could  be  corrected  iterat¬ 
ively)  . 

Finally,  it  may  be  remarked,  that  a  nonlinear  theoretical 
treatment  of  the  problem  would  reveal  more  about  the  low  speed 
oscillations  in  the  resistance  curve  predicted  by  the  1  inoar 
theory . 

Following  Wehausen  (1963),  one  could  assume  the  poten¬ 
tial  to  be  a  power  series  in  terms  of  a  perturbation  para¬ 
meter  such  as  po/oga  or  po/pc2  .  It  should  then  be  fair¬ 
ly  straightforward  to  construct  a  linear  scheme  of  succe  ;s- 
ively  higher  approximations  as  has  already  been  done  for  tl:- 
case  of  a  ship  in  steady  motion.  In  fact,  there  is  rood 
reason  to  believe  that  due  to  the  absence  of  the  hul1  bound¬ 
ary  condition,  the  numerical  evaluation  of  the  resulting  in¬ 
tegrals  here  would  be  somewhat  simpler  than  for  a  ship. 
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Fig.  22  Shallow  Water  Unsteady  Free  Surface  Elevation  (2D) 
for  Fig.  21,  (a)  For  t/g/a  =  0(2)18 
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Fig.  26  Velocity  Pattern  in  Deep  Water  (Constant  Thrust) 
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Fig.  29  Drag  Components  in  Deep  Water  {Constant.  Prop.  Revs.) 
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Fig.  32  Velocity  Pattern  in  Finite  Depth  (Constant  Propeller  Revolutions) 
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This  report  is  cmccrt-ed  vith  the  theoretical  vav<  resistance  cf  an 
air-oassica  vehicle  CACVJ  traveling  over  vater  of  urifcm  finite  tr  ;r- 
f ini te  depth,  in  steady  or  unsteady  notion.  Referring  first  to  steady 
notion,  it  is  shown  that  the  unrealistic  oscillations  in  the  vave  resis¬ 
tance  carve  at  lev  Froude  '-criers  fousd  bv  previous  ver/.err  car.  be  eiir- 
mated  by  using  a  smoothed  cot  pressure  distribution  rather  than  one 
vith  sharp  edges  studied  exclusively  in  the  past.  The  vain  result  cf 
unsteady  motion  calculations  is  that  the  pea>  vave  resistance  in  shal¬ 
low  vater,  ever,  in  moderately  accelerated  notion.  is  appreciably  less 
than  the  corresponding  steady-state  value.  In  fact,  cases  have  been 
found  vhere  an  AC*.’  starting  frees  rest  under  the  action  of  a  constant 
thrust  vouid  sees  to  be  ur.ible  to  cross  the  critical  depth  Froude  number 
on  the  basis  of  quasi-steady  estimates  cf  vave  resistance,  while  the 
more  elaborate  unsteady  calculations  shev  that  it  has  sufficient  rover 
to  reach  its  final  supercritical  cruising  speed.  An  interesting  fea¬ 
ture  of  unsteady  motion  is  that  besides  vave  resistance  there  is  an¬ 
other  mechanism  transferring  energy  to  the  free  surface  vhich  is  here 
called  the  dynamic  sustention  power.  Contrary  to  ir.c-iticr.,  the  wave 
resistance  in  unsteady  motion  ever  finite  depth  sometimes  becomes  nega¬ 
tive  at  supercritical  "rc-ude  numbers  before  finally  approaching  zero 
at  infinite  speed. 
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